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Center:  "Application  of  Thermal  Analysis  to  the  Evaluation  of  LOVA 
Characteristics  of  Propellants"] 

[Text]  Summary 

The  aim  of  this  study  is  to  examine  the  applicability  of  thermal  analysis  to 
characterize  ignitability  of  rocket  propellants  and  their  ingredients.  Double¬ 
base  propellants,  consisting  mainly  of  nitrocellulose  (NC)  and  nitroglycerine 
(NG) ,  were  used  as  a  sample  because  there  is  no  general  agreement  in  reported 
thermograms  of  double-base  propellants  by  others.  We  have  revealed  the 
occurrence  of  two  distinctively  different  ignitions  triggered  by  exothermic 
decomposition  of  NC  and  NG:  as  heating  rate  was  increased,  transition  occurs 
from  NC-triggered  ignition  to  NG-triggered  ignition.  We  have  also  demonstrated 
that  both  heating  rate  and  sample  weight  predominate  the  onset  of  ignition;  at 
a  constant  sample  weight,  ignition  occurs  at  a  higher  temperature  as  heating 
rate  is  raised  along  with  the  NG  or  NC  decomposition  peaks;  at  a  constant 
heating  rate,  transition  occurs  from  NC-triggered  ignition  to  NG-triggered 
ignition  as  sample  weight  is  added.  We  also  confirmed  that  the  ignition  of  the 
double-base  propellants  originates  in  the  condensed-phase  under  the  experi¬ 
mental  conditions  tested.  The  thermonanalytical  method  has  been  proven  to  be 
a  promising  diagnostic  technique  to  characterize  heat-resistant  properties  of 
propellants  and  explosives  by  virtue  of  simplicity  of  data  acquisition  and 
analysis . 

1.  Introduction 

The  importance  of  creating  rocket  propellants  and  explosives  that  are  not 
sensitive  to  heat  and  impact  has  been  a  common  subject  in  the  field  of  guided 
weapons,  firearms,  and  ammunition.  The  double— base  propellant  (main 
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components:  nitrocellulose  and  nitroglycerine)  is  commonly  used  as  a  smokeless 
rocket  propellant  or  as  a  gun  propellant  powder,  and  it  appears  that  the 
performance  of  the  double-base  propellant  can  be  enhanced  without  degrading 
its  characteristic  smokelessness.  Therefore,  in  order  to  improve  the  heat 
resistance  of  this  kind  of  propellant,  i.e.,  to  develop  low  vulnerability 
(LOVA)  propellants,  it  is  absolutely  necessary  to  obtain  basic  data  on  the 
ignition  characteristics  and  behavior  of  double-base  propellants  when 
subjected  to  heat. 

The  purpose  of  this  research  work  was  to  use  a  thermal  analyzer  to  evaluate 
the  ignition  characteristics  and  to  clarify  the  relationship  between  the 
ignition  characteristics  and  the  exothermic  physicochemical  characteristics 
that  accompany  thermal  decomposition.  Also,  this  thermal  analyzer  is  used  to 
measure  changes  in  the  mass  of  the  samples  and  the  input  and  output  of  heat 
caused  by  raising  the  temperature  of  these  samples  at  a  constant  heating  rate. 
Thermal  analyzer  has  been  applied  to  the  evaluation  of  the  safety  of  hazardous 
materials,  and  there  are  some  research  reports  on  this  application,  but  there 
is  no  example  of  a  study  where  ignition  characteristics  have  been  subjected  to 
a  quantitative  evaluation.  We  have  studied  the  possibility  of  using  a  thermal 
analyzer  to  clarify  ignition  mechanisms  and  to  evaluate  ignition  characteris¬ 
tics  ,  thereby  taking  advantage  of  the  fact  that  it  is  easy  for  a  thermal 
analyzer  to  change  variables  that  affect  ignition,  such  as  ambient  gas, 
pressure,  etc.,  in  sample  rooms. 

Many  researchers  have  long  been  engaged  in  research  on  the  thermal  decomposi¬ 
tion  characteristics  of  the  primary  components  of  double-base  propellants; 
nevertheless,  they  continue  to  disagree  with  each  other  on  some  important 
points.^  In  our  research  laboratory,  we  have  carried  out  research  on  double¬ 
base  propellants  and  the  thermal  decomposition  speed  and  decomposition 
mechanism  of  nitrates,  which  are  the  principal  components  of  these  double-base 
propellants.  Our  goal  has  been  to  resolve  a  number  of  unsolved  problems.  We 
have  determined  that  an  autocatalyzed  reaction  is  generated  at  a  rate¬ 
determining  step  of  thermal  decomposition  in  comparatively  high  temperature 
regions  and  that  this  autocatalyzed  reaction  is  an  oxidative  reaction  produced 
by  the  nitric  oxide  generated  by  thermal  decomposition.^’^  This  report 
describes  the  results  of  our  efforts  to  develop  simple  methoas  for  evaluating 
the  thermal  ignition  characteristics  necessary  for  improving  LOVA  character¬ 
istics  against  heat  in  light  of  the  results  of  the  basic  research  work  that 
has  been  carried  out  up  to  now. 

2 .  Experiment 

2 . 1  Sample 

Three  kinds  of  double-base  propellants  were  used  in  this  experiment;  the  ratio 
of  their  main  components,  i.e.,  nitrocellulose  (NC)  to  nitroglycerine  (NG) 
range  from  about  1.4: 1-5:1.  The  samples  were  formed  into  cylinders  with  a 
diameter  of  1  mm  and  a  height  of  1  mm.  Table  1  shows  the  composition  of  the 
samples.  These  samples  are  all  so-called,  "Non-Plateau  Double-Base 
Propellant,"  which  means  they  do  not  contain  any  combustion  catalysts  such  as 
lead  salt,  etc. 
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Table  1.  Composition  of  Propellant  Samples 


Components 

Propellant  1 

Propellant  2 

Propellant  3 

Nitrocellulose  (NC) 

51.8 

71.2 

79.3 

Nitroglycerin  (NG) 

36.5 

24.1 

16.0 

DEP  (Diethyl  phthalate) 

9.6 

2.4 

2.4 

2-NDPA  (2-nitro  diphenyl 
amine ) 

2.1 

2.1 

2.1 

NC/NG 

1.4 

3.0 

5.0 

H2C-O-NO2 

HC-O-NO2 

H2C-O-NO2 


Nitroglycerin 


2.2  Devices 

A  differential  scanning  calorimeter  (DSC— 20B)  and  a  thermal  balance  (TGA— 20B) 
with  an  analog  differentiator  were  used  as  thermal  analyzers.  A  sample  was  put 
in  an  aluminum  sample— receiving  dish  with  a  diameter  of  5.8  mm  and  a  height  of 
1.8  mm  in  a  releasing  state,  and  the  influence  was  investigated  while  changing 
the  weight  of  the  sample  across  a  range  of  about  0.2-10  mg.  The  thermal 
balance  is  of  the  vertical  hanging  type  and  is  structured  so  that  a  sample 
receiving  dish  is  put  under  a  platinum  hanging  wire.  Temperatures  were 
calibrated  by  measuring  the  magnetic  transition  point  (Curie  point:  353°C)  of 
nickel . 

Figure  1  shows  the  structure  and  principle  of  operation  of  the  thermal  balance 
(TGA-20B)  and  of  the  differential  scanning  calorimeter  (DSC-20B) .  Both  devices 
are  DT-20  series  devices  made  by  Shimadzu  Corp.  and  are  so-called  "Micro- 
Thermal  Analyzers."''  They  can  be  used  to  measure  the  weight  of  samples  as  small 
as  several  milligrams.  A  preliminary  test  was  conducted  by  using  a  mixture  of 
nitrogen  and  air  as  the  atmosphere  for  the  measurement.  As  a  result,  it  was 
ascertained  that  the  thermal  ignition  of  double-base  propellants  is  not 
affected  by  any  partial  pressure  of  oxygen  (Figure  14).  Accordingly,  the 
measurement  was  carried  out  in  untreated  air  without  any  use  of  inactive  gas . 
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Simultaneous  TG-DTG  Apparatus 
(Shimadzu  TGA-20) 


Figure  1.  Block  Diagrams  of 
3.  Result  and  Discussion 


BELL  JAR 


Differential  Scanning  Calorimeter 
(Shimadzu  DSC'‘208) 


and  TG— DTG^  Used  in  This  Research 


3.1  Relationship  Between  Thermal  Decomposition  Characteristics  and  Thermal 
Firing  Properties 

There  have  been  many  reports  on  thermal  analyses  of  double-base  propellants. 
However,  the  decomposition  peaks  of  NG  and  NC,  which  are  the  principal 
components  of  double-base  propellants,  have  not  been  separated  from  each  other 
by  thermal  measurement  using  differential  thermal  analysis  (DTA)  or  a 
f fstential  scanning  calorimeter  (DSC) ,  nor  have  they  undergone  weight 
reducing  measurement  based  on  thermal  weight  analysis  (TG).^“®  It  is  believed 
that  this  is  due  to  the  following  factors:  1)  both  thermal  decomposition 
temperatures  are  very  close  to  each  other;  2)  both  thermal  decomposition  speed 
characteristics  are  similar  to  each  other;  and  3)  the  decomposition  peaks 
cannot  be  separated  because  the  thermal  analyzer  does  not  have  sufficient 
thermal  response.  However,  it  is  difficult  to  believe  that  the  third  factor  is 
the  principal  problem,  because  the  two  peaks  cannot  be  separated  from  each 
other  even  by  the  DSC,  which  has  a  very  high  thermal  response.  Another 
possible  cause  is  that  there  exists  a  number  of  small  approaching  exothermic 
peaks  that  are  hidden  because  another  heating  value  is  large. 

In  this  research  work,  then,  a  weight  loss  speed  curve  was  found  by  using  a 
highly  sensitive  Kern-type  (vertical  system)  thermo-balance.  This  was  compared 
with  a  heat  generation  speed  curve  determined  by  using  the  DSC.  In  addition, 
as  shown  in  Table  1,  samples  with  NC/NG  ratios  of  about  3:1  and  5:1,  as  well 
as  a  standard  propellant  with  a  ratio  of  almost  1:1  were  used  to  ensure  the 
separation.  As  a  result,  in  a  typical  example,  as  shown  in  Figure  2,  it  was 
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possible  to  clearly  distinguish  two  peaks  by  a  differential  thermal  weight 
analysis  (DIG)  curve.  Reproducibility  is  low,  but  an  exothermic  shoulder  may 
be  recognized  at  the  low  temperature  side  in  the  DSC  curve. 


Figure  2.  Result  of  Thermochemical  Analysis  of 
Double— Base  Propellant  (TG— DTG:  DSC) 

Further,  because  samples  with  different  NC/NG  ratios  were  used,  it  was 
possible  to  identify  the  first  peak  and  the  second  peak.  That  is,  it  has  been 
determined  that  the  first  peak  is  due  to  the  decomposition  of  NG  while  the 
second  peak  is  due  to  that  of  NC.  This  is  because  the  height  of  the  first  peak 
is  inversely  proportional  to  the  NC/NG  ratio  (Figure  3) . 


Figure  3.  Comparison  of  TG-DTG  Curves  for  Double-Base 
Propellants  With  Different  NC/NG  Ratios 
This  figure  shows  that  the  first  peak  (low  temperature  side)  and 
the  second  peak  are  caused  by  the  thermal  decomposition  of  NG  and 
of  NC,  respectively. 
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The  relationship  between  a  propellant's  thermal  decomposition  characteristics 
and  its  thermal  ignition  properties  was  clarified  by  changing  the  weight  of 
the  samples  and  the  rate  at  which  they  were  heated.  It  is  possible  to  clearly 
distinguish  the  occurrence  of  ignition,  because  the  weight  apparently 
increases  momentarily  at  the  time  of  ignition,  i.e. ,  this  increase  appears  as 
an  upward  spike  on  the  DTG.  As  shown  in  Figure  4,  ignition  is  generated 
qualitatively  in  the  vicinity  of  temperature  close  to  the  two  thermal 
decomposition  peaks,  and  when  the  weight  increases  at  a  constant  heating  rate, 
the  ignition  will  transit  from  the  first  peak  to  the  second  peak.  This 
transition  phenomenon  was  not  previously  known. 


SO  100  150  200  250 
TEMPEfiATURE,  ‘C 

Figure  4.  TG-DTG  Curve  for  Double— Base  Propellant 

When  the  Weight  of  the  Samples  Is  Changed 

Figure  5  shows  the  correlation  between  the  rate  of  weight  reduction  at  firing 
time  and  ignition  temperature  for  three  kinds  of  samples.  The  ignition 
temperature  does  not  depend  on  the  NC/NG  ratio,  but  the  rate  of  weight 
reduction  at  the  time  of  ignition  greatly  depends  on  the  NC/NG  ratio.  It  was 
learned  that  the  ignition  transited  at  a  weight  reducing  rate  equivalent  to 
the  sum  of  the  plasticizer  (DEP)  and  the  NG,  which  is  an  almost  liquid 
component,  present  in  the  propellant. 

In  this  way,  when  the  number  of  samples  was  small  at  a  constant  heating  rate, 
only  decomposition  occurred,  while  thermal  ignition  at  the  second  stage  was 
separated  and  measured  together  with  increase  in  the  weight  of  the  samples. 
It  seems  that  this  can  be  explained  by  a  fortuitous  set  of  experimental 
conditions.  Next,  we  will  discuss  the  conditions  necessary  for  two— stage 
ignition  and  the  influence  of  heating  rate  on  thermal  ignition. 

The  weight  of  typical  double-base  propellants  (characterized  by  propellant  1) 
as  samples  is  selected  properly  under  atmospheric  pressure.  Figure  6  shows  the 
result  of  the  selection.  As  shown  in  this  figure,  thermal  decomposition  is 
represented  by  two  straight  lines  such  that  when  the  logarithm  of  the  heating 
rate  is  taken  on  the  x-axis  and  the  reciprocal  number  of  the  decomposition 
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Figure  5.  Influence  of  NC/NG  Ratio  on  Second-Stage  Ignition 

This  figure  shows  that  the  first  ignition  peak  is  generated  at  a 
value  close  to  the  total  sum  of  the  content  of  NG  and  DEP,  and  is  caused 
by  the  exothermic  decomposition  of  NG. 


PEAK  TEMERATURE  OF  MASS-LOSS  PJ\TE(DTG),  °C 
200  190  180  170 


Figure  6.  Characteristics  of  Thermal  Decomposition  Ratios  for  NG  and  NG, 

the  Primary  Components  of  Double— Base  Propellants ,  and  Influence 
of  Pressure  on  NC  and  NG 
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peak  temperature  is  taken  on  the  transverse  axis,  the  gradients  (activation 
energy)  of  the  two  straight  lines  will  differ  from  each  other.®  The  straight 
line  on  the  low  temperature  side  is  due  to  the  thermal  decomposition  of  NG, 
and  the  apparent  activation  energy  is  abnormally  high,  reaching  about  60 
Kcal/mol.  Meanwhile,  the  straight  line  at  the  high  temperature  side  is  due  to 
the  thermal  decomposition  of  NC,  and  the  activation  energy  is  about 
47  Kcal/mol.  Variations  can  be  seen  in  ignition  temperatures,  but  thermal 
ignition  is  arranged  on  a  straight  line  that  almost  equals  thermal 
decomposition  and  transits  from  the  second  peak  ignition  (NC  firing)  to  the 
first  peak  ignition  (NC  firing)  in  accordance  with  an  increase  in  the  heating 
rate  (Figure  7) . 


PEAK  TEMPERATURE  OF  MASS-LOSS  PJ^TE,  ^0 


Figure  7 .  Relation  Between  Thermal  Decomposition  Characteristics  and 

Transition  Phenomena  of  the  Second  Stage  Ignition  of  Double- 
Base  Propellant 

Many  values  for  the  activation  energies  of  the  thermal  decomposition  of 
double-base  propellants  have  been  reported.  These  can  be  broadly  classified 
into  two  values,  i.e.,  about  35  Kcal/mol  and  about  45-50  Kcal/mol.  It  has  been 
considered  up  to  now  that  the  former  value  is  due  to  an  autocatalyzed  reaction 
while  the  latter  value  is  due  to  homo  lysis.  ^  In  our  research  laboratory, 
however,  we  recently  reached  a  conclusion  quite  the  reverse  of  the 
conventionally  established  theory.  This  is  based  on  the  results  of  our 
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research  on  the  thermal  decomposition  of  double-base  propellants  at 
comparatively  low  temperatures  (150°C  or  less) ,  and  we  have  reported  on  this 
conclusion.  In  other  words,  a  logical  analysis  of  the  reaction  rate  has 
demonstrated  that  a  value  of  about  35  Kcal/mol  is  due  to  the  homolysis  of  the 
combination  of  N-0,  while  the  activation  energy  close  to  50  Kcal/mol  is  due  to 
the  autocatalyzed  reaction. It  was  also  anticipated  that  the  main  reaction 
(rate-determining  step)  in  high  temperature  regions  would  be  autocatalyzed 
reaction  with  high  activation  energies .  As  shown  in  Figure  7 ,  a  value  of 
47  Kcal/mol  is  the  activation  energy  for  the  thermal  decomposition  of  NC 
contained  in  a  double-base  propellant  in  air  at  1  atmosphere  of  pressure. 
This  value  is  consistent  with  the  expected  value. 

Meanwhile,  the  value  of  the  activation  energy  for  the  decomposition  of  NG  at 
1  atmosphere  of  pressure  was  abnormally  high,  reaching  about  60  Kcal/mol.  As 
shown  in  Figure  6,  that  of  NG  at  low  pressure  (0.16  atmosphere)  was  as  low  as 

46  Kcal/mol.  Likewise,  in  the  same  way  as  at  1  atmosphere,  that  of  NC  was 

47  Kcal/mol.  It  was  unknown  before  now  that  the  activation  energy  for  the 
thermal  decomposition  of  nitrate  is  dependent  on  pressure  to  such  an  extent. 
It  is  not  clear  why  only  the  NG  contained  in  double-base  propellants  is 
affected  by  pressure.  This  must  be  determined  by  future  research. 

The  introduction  to  this  section  indicates  that  there  has  been  no  reporting  up 
to  now  on  the  fact  that  the  thermal  analyzer  curve  of  double-base  propellants 
was  divided  into  two  peaks.  It  was  believed  that  this  was  because  one  heating 
value  was  larger  than  the  other  and  concealed  the  second.  In  order  to  confirm 
this  hypothesis,  the  heating  values  of  NG  and  NC  were  determined  based  on  a 
comparison  of  the  weight  reducing  speed  (DSC)  curves  obtained  from  this 
research  work  on  the  assumption  that  the  heating  value  is  a  proportional 
constant  of  the  heat  generation  speed  and  weight  reducing  speed.  That  is,  as 
shown  in  equation  (1),  the  heating  value  (AHr)  at  an  arbitrary  temperature  is 
a  proportional  constant  of  the  height  of  the  DSC  curve  and  that  of  the  DTG 
curve . 


AHr=  (dj,)  T-l  (d„)  r 

Height  of  DSC  curve  {meal/ sec) 

Height  of  DTG  curve  {mg/ sec) 

The  advantage  of  this  equation  (differential  equation)  is  that  the  heating 
value  can  be  found  as  a  function  of  temperatures .  Figure  8  shows  examples  of 
the  results  of  the  calculations  and  the  DTG  and  DSC  curves  of  a  standard 
double-base  propellants  (propellant  1).  The  heating  value  of  NG  and  NC  were 
found  to  be  about  50-60  and  about  200  cal/g,  while  that  of  double-base 
propellants  is  about  70-80  cal/g. These  values  are  almost  identical  to  those 
measured  in  this  research  work. 

The  measured  value  for  the  decomposition  heating  value  of  NC  is  equivalent  to 
20  percent  of  its  explosive  heating  value  (about  1,000  cal/g), while  that  of 
NG  is  equivalent  only  to  about  3  percent  of  its  explosive  heating  value 
(1,500  cal/g). “ 
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Figure  8.  Thermal  Dependence  of  Value  of  Heat  Generated  by  Thermal 
Decomposition  of  Double— Base  Propellant 

The  values  of  heat  generated  from  NG  to  NG  were  measured  to  be 
50—60  cal/g  and  about  200  cal/g,  respectively. 


The  relationship  between  explosive  heating  value  and  heating  value  is  not 
proportional.  It  is  believed  that  this  fact  may  provide  some  clues  about  the 
decomposition  mechanism.  The  value  obtained  by  measuring  the  decomposition 
heating  value  of  NG  includes  errors,  because  an  endothermic  reaction  is  caused 
by  evaporation  at  the  same  time  as  that  of  exothermic  decomposition. 
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Therefore,  it  is  necessary  to  pay  attention  to  the  fact  that  the  measured 
value  is  lower  than  the  actual  value . 

As  a  result  of  the  study  described  in  section  3.3,  it  was  determined  that  the 
exothermic  influence  in  gaseous  phases  can  be  ignored.  Therefore,  it  can  be 
considered  that  the  values  obtained  by  measuring  these  decomposition  heating 
values  are  almost  purely  heating  values  in  condensation  phases. 

We  will  now  move  on  to  discuss  the  difference  between  the  heating  value  of  NC 
and  that  of  NG.  The  difference  between  the  two  chemical  structures  is  as 
follows:  the  former  is  a  high  molecular  compound  with  a  repeating  unit  in  a 
pyranose  ring  (grape  sugar  structure)  where  oxygen  atoms  are  contained  in  six 
member  rings,  while  the  latter  is  a  nitrate  with  a  low  molecular  weight  where 
glycerin,  which  is  a  trivalent  alcohol,  is  nitrated. 

In  the  case  of  NG,  it  is  believed  that  aldehydes  (for  example,  formaldehyde: 
HCHO)  with  a  low  molecular  weight  are  produced  and  are  separated  readily  from 
reaction  systems  following  the  release  of  nitrogen  dioxide  due  to  homolysis  of 
the  weakest  bonding  of  oxygen  and  nitrogen.  Therefore,  it  is  felt  that 
nitrogen  dioxide  is  not  a  sufficient  cause  for  any  oxide  reaction  at  low 
pressures,  such  as  atmospheric  pressure.  We  will  quantitatively  discuss  this 
point  in  the  next  item  and,  as  a  conclusion,  we  have  determined  that  oxidation 
due  to  NO2  cannot  occur  within  the  conditions  of  this  experiment.  On  the  other 
hand,  in  the  case  of  NC,  it  is  known  that  homolysis  of  nitrate  group  (-O-NO2) 
side  chains,  homolysis  of  the  pyranose  ring,  and  cleavage  of  the  principal 
chain  occur  simultaneously  and  that  complex  metathetical  reactions  occur  in 
the  viscous  condensation  phases  until  the  NC  turns  to  a  fragment  with  a 
molecular  weight  sufficient  for  evaporation.^  Consequently,  it  can  be  said  that 
oxidation  is  induced  by  NO2  and  increases  in  the  heating  value.  Actually,  the 
results  we  obtained  indicate  that  the  thermal  decomposition  peak  temperature 
of  NC  is  about  10°C  higher  than  that  of  NG  (Figures  3  and  6).  It  has  also  been 
demonstrated  that  the  thermal  decomposition  peak  temperature  of  NC  decreases 
along  with  a  decrease  in  pressure.  Thus,  as  shown  in  Figure  6,  the  thermal 
decomposition  peak  temperature  of  NC  is  about  3-4°C  lower  than  the  peak  value 
at  0.16  atmosphere. 

It  was  found  in  this  research  work  that  the  double-base  propellant  can  cause 
a  two— stage  ignition  depending  on  the  heating  rate  and  the  number  of  samples. 
According  to  one  thermal  ignition  theory,  this  phenomenon  can  be  explained  as 
follows : 

It  is  believed  that  when  the  pace  of  heat  generation  surpasses  that  of  heat 
loss,  thermal  ignition  will  occur.  Now,  considering  the  temperature  differ¬ 
ences  and  the  temperatures  on  the  ordinate  and  transverse  axes,  respectively, 
as  shown  in  Figure  9,  in  general,  the  rate  of  heat  loss  can  be  expressed  by  a 
straight  line  that  passes  through  the  origin.  The  inclination  of  this  straight 
line  is  peculiar  to  thermal  analyzers.  In  the  case  of  double-base  propellants, 
the  rate  of  heat  generation  can  be  expressed  with  a  curve  that  has  two  peaks 
(see  the  DSC  curve  shown  in  Figure  2).  With  respect  to  the  conditions  govern¬ 
ing  the  occurrence  of  two— stage  ignition,  the  following  three  cases  can  be 
considered  depending  on  whether  or  not  the  straight  line  representing  heat 
loss  and  the  heat  generating  curve  intersect. 
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TEMPERATURE  - ► 


Figure  9.  Characteristic  Representation  of  How  Conditions  Governing  Second 
Stage  Ignition  Are  Explained  by  Theories  of  Thermal  Ignition 
As  can  be  seen  from  this  figure,  second  stage  ignition  occurs  when 
the  first  peak  (low  temperature  side)  is  lower  than  the  second  (high 
temperature  side) ,  and  the  ignition  temperature  transits  from  the  second 
peak  to  the  first  one.  The  result  predicted  by  this  figure  accords  with 
the  experimental  results  we  obtained. 

In  the  first  case  shown  in  Figure  9,  thermal  ignition  does  not  occur,  because 
the  heat  loss  straight  line  is  higher  than  the  heat  generation  curve.  In  the 
second  case  shown  in  the  figure,  there  is  a  possibility  that  thermal  ignition 
could  occur  at  temperatures  higher  than  the  temperature  of  the  intersection 
point,  because  the  heat  loss  straight  line  intersects  the  second  peak  of  the 
heat  generation  curve.  Finally,  in  the  third  case  shown  in  the  figure,  there 
is  a  possibility  that  thermal  ignition  could  occur  in  the  vicinity  of  the 
first  peak,  because  the  heat  loss  straight  line  intersects  the  first  peak. 

Assuming  that  the  former  represents  first  peak  ignition  and  the  latter  second 
peak  ignition,  the  first  and  second  peak  ignitions  are  based  on  the  exothermic 
decomposition  of  NG  and  NC,  respectively.  A  comparison  of  cases  2  and  3 
clearly  indicates  a  condition  for  generating  two-stage  ignition.  That  is,  the 
first  exothermic  peak  on  the  low  temperature  side  is  lower  than  the  second 
peak  on  the  high  temperature  side. 

3.2  Influence  of  Weight  of  Samples 

The  previous  section  outlined  the  influence  of  the  weight  of  samples  in  the 
ignition  characteristics  of  double-base  propellants  (Figure  4) .  This 
Influence  will  now  be  studied  in  greater  detail.  Figure  10  shows  the  result  of 
an  investigation  of  the  effect  of  the  weight  reducing  rate  at  ignition  by 
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changing  the  weight  of  sample  2  (NG  content  =  about  24.1  percent)  from  about 
2-8  mg.  Second  peak  ignition  was  generated  in  a  range  of  sample  weights  up  to 
about  4  mg,  while  first  peak  ignition  was  generated  at  weights  of  about  4  mg 
or  more.  The  weight  loss  is  about  24  percent  during  a  transition  from  second 
peak  ignition  to  first  peak  ignition.  This  values  accords  with  NG  content. 
This  experimental  result  has  demonstrated  that  first  peak  ignition  is  caused 
by  NG. 


Figure  10.  Influence  of  Sample 
Weights  on  Second- 
Stage  Ignition  (Constant 
heating  rate) 


Figure  11.  Influence  of  Weight  of 

Samples  and  Heating  Rate 
on  Ignitability  Limit 


Figure  11  shows  the  result  of  plotting  logarithms  of  the  weights  of  samples 
and  heating  rate  in  order  to  clarify  the  extent  of  the  influence  of  the  weight 
of  these  samples  and  the  heating  rate  on  two— stage  ignition.  It  shows  that 
ignition— generating  critical  areas  have  an  inclination  of  about  minus  one- 
half,  that  the  weight  of  the  samples  and  the  heating  rate  affect  ignition 
almost  simultaneously,  and  that  weight  and  heating  rate  are  in  inverse 
proportion  to  each  other.  In  other  words,  when  the  thermal  ignition 
characteristics  of  explosives  are  assessed  by  using  a  thermal  analyzer,  it  is 
necessary  to  control  the  weight  of  the  samples  and  the  heating  rate  exactly. 


3.3  Influence  of  Pressure  and  Oxygen 

Earlier  studies  have  shown  that  the  thermal  ignition  of  double-base 
propellants  is  caused  by  exothermic  reactions  of  the  main  components  in 
condensation  phases.  An  experiment  on  the  influence  of  pressure  was  conducted 
to  verify  that  ignition  is  in  fact  generated  during  such  condensation  phases. 
This  section  will  describe  the  design  and  results  of  this  experiment. 


In  order  to  study  the  influence  of  pressure  theoretically,  it  is  necessary  to 
understand  the  reaction  mechanism  for  the  thermal  decomposition  of  nitrate. 
However,  many  aspects  of  this  reaction  mechanism  are  unclear,  and  we  are 
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currently  working  to  clarify  the  reaction  mechanism  in  our  research 
laboratory.  So  far,  we  have  demonstrated  by  using  a  chemiluminescence  method 
and  an  infrared  spectral  method  that  the  above  reaction  mechanism  is  an  oxygen 
reaction  in  which  peroxides  and  hydrogen  peroxides  are  intermediate  products . 
We  will  not  attempt  to  describe  the  decomposition  mechanism  in  detail,  but 
will  carry  out  macroscopic  studies  based  on  the  measurement  of  the  heating 
values  of  NG  and  NC  described  in  the  previous  section. 

There  is  a  substantial  body  of  research  on  the  combustion  of  double-base 
propellants,  and  it  has  been  learned  that  the  structure  of  the  combustion 
waves  and  main  reactions  are  like  those  shown  in  Figure  12.  Reaction  1  is  a 
homolysis  of  the  bonding  of  oxygen  and  nitrogen  in  a  nitrate  (RONO2) ;  it  occurs 
in  the  condensation  phases.  Reaction  2  is  an  exothermic  reaction  of  nitrogen 
dioxide  (NO2)  and  aldehydes;  it  occurs  in  the  condensation  phases  or  in  phases 
very  close  to  these  condensation  phases.  Reaction  3  is  an  oxygen  reaction  in 
which  flames  are  generated  by  the  produced  nitric  oxide;  it  does  not  occur  in 
any  inactive  gases  at  atmosphere  pressure  and  is  called  "Fitz  Burning."  The 
heating  value  in  the  case  of  reaction  3  is  almost  half  (about  500  cal/g)  of 
the  explosion  heating  value.  Therefore,  it  is  believed  that  all  that  is 
necessary  is  to  consider  reaction  2  in  experimental  conditions  at  pressures 
lower  than  atmospheric  pressure  and  at  atmospheric  pressures  in  nitrogen 
atmospheres  used  in  this  research  work. 


Reaction  1  Reaction  2  Reaction  3 

NO,  - ^  NO.  H*  - ►  Nz,  CO 

R’ClIO  CO,  CO,  COi.  HzO 

Nonluminou8  Luminous 

Fizz  phase  flame  phase  flame  phase 

'  (Gaseous 

phase ) 

Figure  12.  Principal  Reactions^  and  Structure  of 

Combustion  Waves  of  Double-Base  Propellant 

The  chemical  structural  formula  of  NG  is  represented  by  (CH20N02)2  (CHONO2) . 
Now,  let  us  suppose  that  3  mols  of  formaldehyde  (HCHO)  and  3  mols  of  nitrogen 
dioxide  are  produced  from  1  mol  of  NG,  respectively,  in  order  to  simplify  this 
formula.  There  has  been  considerable  research  on  the  vapor  phase  reactions  of 
formaldehyde  and  nitrogen  dioxide,  and  it  is  known  that  reaction  rate 
parameters  (activation  energy  and  pre-exponential  factor)  change  at 
temperatures  around  160 °C  in  the  secondary  reaction. 

The  following  provides  a  comparison  of  the  diffusion  rate  of  atmospheric  gas 
and  the  reaction  rate  under  the  conditions  of  this  experiment  carried  out  by 
taking  the  characteristic  time  (rcj,)  as  a  parameter.  Assximing  that  the 
characteristic  reaction  time  is  (rj-j,) .  the  half-life  is  (.^1/2} ,  the  reaction 
rate  constant  is  k  (1  mol~^second~^)  ,  and  the  initial  concentration  is  a 
(mol  1  ^) .  From  this,  the  following  equation  can  be  obtained: 


RONOz  - ► 

r  0  p  e  1 1  a  n  t;:;;;:;:;:;:; 

•X;;:;(Condensa  - 
•^Vtion  phasej-:;:;:;:;; 
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T^C*=T^l/2  =  -L/-K-a 


u; 


Now,  when  k  =  10^°e  i9.ooo/RT  (-p  >  160°C)  ,  it  is  estimated  that  a  value  of 
k  =  1.2  (1  mol“^  second  “^)  can  be  obtained  at  a  temperature  of  200°C  and  that 
a  value  of  a  =  10“^ (mol  1”^)  can  be  obtained  at  atmospheric  pressures.  There¬ 
fore,  Teh  is  about  10  seconds.  Meanwhile,  the  characteristic  time  (tq)  for 
molecular  diffusion  can  be  represented  by  the  equation; 

To=dV8D  (3) 

Where,  d  (centimeter)  =  width  of  samples,  and 

D  (square  centimeters  per  second)  =  a  diffusion  constant 

a  decrease  in  pressure  will  cause  a  decrease  in  a  and  an  Increase  in  the  value 
of  Teh  in  equation  (2).  However,  a  decrease  in  pressure  will  increase  the 
value  of  D  and  will  decrease  the  value  of  ro  in  equation  (3).  Therefore,  the 
value  of  Tch/ro  is  proportional  to  the  second  power  of  pressure,  and  with 
further  decreases  of  this  pressure,  the  possibility  of  ignition  of  gaseous 
phases  is  further  reduced. 

Now,  substituting  the  diameter  (d  =  0.1  cm)  of  the  samples  used  in  this 
experiment  and  the  diffusion  constant  (D  =  0.5  cmVs)  for  equation  (3),  we  have 
Tp  =  10~^  seconds.  Therefore,  the  characteristic  time  of  vapor  phase  reaction 
at  atmospheric  pressures  is  about  10''  times  longer  than  that  of  diffusion,  the 
exothermic  oxidizing  reaction  (reaction  2  shown  in  Figure  13)  can  be  ignored 
in  gaseous  phases  of  aldehyde  due  to  the  NO2,  and  the  possibility  of  ignition 
in  the  gaseous  phases  can  be  denied. 


Figure  13 .  Influence  of  Pressure  on  Ignition 
The  absence  of  pressure  suggests 
ignition  of  the  condensation  phase. 


For  example,  when  the  diameter  of  the  samples  is  10  times  that  of  the  above 
and  the  pressure  is  150  atmospheres  both  characteristic  times  for  the  systems 
used  in  this  experiment  will  be  equalized.  From  the  results  obtained  by 
comparing  and  studying  the  characteristic  times,  we  concluded  that  ignition 
does  not  depend  on  any  pressure  and  partial  pressure  of  oxygen.  This  conclu¬ 
sion  is  supported  by  the  measurement  results  we  obtained  (Figures  13  and  14)  . 
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Figure  14.  Influence  of  Oxygen  Concentration 
on  Ignition 

Ignition  is  not  affected  by  any  amount  of  oxygen; 
this  suggests  ignition  at  the  condensation  phase. 
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4.  Conclusion 

We  have  carried  out  research  on  the  possibility  of  using  thermal  analyzers  to 
evaluate  LOVA  characteristics  during  low— speed  heating  together  with  the 
thermal  ignition  characteristics  of  explosives.  We  have  confirmed  that  these 
thermal  analyzers  are  very  effective  for  quantitatively  describing  the  thermal 
decomposition  characteristics  of  explosives.  Double-base  propellants  are 
widely  used  as  rocket  propellants  and  propelling  charges.  We  have  discovered 
the  following  new  facts  by  using  such  double-base  propellants  as  samples  in 
this  research  work. 

1.  The  principal  components,  NG  and  NC,  can  be  distinguished  by  analyzing 
their  thermal  weights  (TG  and  ?);  they  appear  as  exothermic  decomposition 
peaks  adjacent  to  each  other. 

2.  The  ignition  of  double-base  propellants  during  low-speed  heating 
(1-20° C/minute)  occurs  in  the  condensation  phases.  The  ignition  temperature 
depends  on  the  heating  rate  and  the  weight  of  the  samples,  and  is  close  to  the 
exothermic  decomposition  peak  temperature  of  the  main  components,  NG  and  NC. 
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3.  When  the  weight  of  the  samples  is  constant  and  the  heating  rate  is  raised, 
ignition  due  to  NC  will  transit  to  that  due  to  NG,  which  is  about  10°C  lower 
than  the  temperature  of  the  ignition  due  to  NC . 

Considering  the  results  that  have  been  described  above,  and  the  advantage 
whereby  important  thermochemical  parameters  such  as  heating  values ,  etc .  ,  can 
be  ascertained  from  thermal  analyzing  measurements,  it  can  be  stated  as  a 
conclusion  that  the  thermal  analyzer  is  useful  for  evaluating  complex  ignition 
characteristics  and  for  clarifying  ignition  phenomena. 

Also,  it  is  absolutely  necessary  to  clarify  the  conditions  governing  ignition 
(runaway  chemical  reaction)  in  order  to  clarify  pure  thermal  decomposition 
reactions  of  propellants  that  are  easy  to  ignite,  explosives,  and  detonating 
explosives  by  using  a  thermal  analyzer. 
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[Article  by  Takao  Komi,  Masakatsu  Sugii,  and  Hiroshi  Kara,  assistant  directors 
general  of  the  Second  Research  Laboratory  of  the  Third  Optical  Wave  Department 
of  TDRI's  Second  Research  Center:  "TRDI  Technical  Report  on  Ar-Xe  Laser"] 

[Text]  Summary 

This  paper  presents  the  laser  performance  of  Ar-Xe  TE  laser  with  UV- 
preionizers.  We  obtained  a  maximum  output  power  of  1.4  mJ  at  a  gas  composition 
of  Ar:Xe  =  100:1,  total  pressure  of  1  atmosphere  with  a  charged  voltage  of 
25  kV,  a  main  condenser  capacity  of  60  nF,  and  a  peaking  condenser  capacity  of 
54  nF. 

1.  Introduction 

The  Ar-Xe  laser  employs  photons  that  are  discharged  when  electrons  transit 
from  the  5d  orbit  to  the  6p  orbit  of  an  excited  Xe  atom.  These  photons  have  a 
wide  oscillation  wavelength  of  1-3  in  the  near  infrared  region.  As  for 

this  wavelength,  a  wavelength  of  1.73  ^m  is  said  to  have  a  comparatively  high 
intensity.  A  laser  with  this  wavelength^  possesses  the  following  characteris¬ 
tics:  1)  its  ability  to  penetrate  the  atmosphere  is  comparatively  good;  and 
2)  it  is  safe  for  the  eyes.  Also,  this  laser  theoretically  should  have  a  long 
life,  because  it  employs  only  inert  gases. 

Experimental  research  on  such  a  laser  has  been  conducted  since  around  1970  and 
various  excitation  methods^”®  have  been  devised,  but  the  laser  is  liable  to 
experience  arc  discharge  and  it  becomes  difficult  to  obtain  uniform  glow 
discharge,  accompanied  by  high  atmospheric  pressure,  because  the  ionization 
potential  of  the  service  gas  is  comparatively  low.  At  present,  high  output 
and  high  efficiency®  can  be  obtained  from  electric  discharge  using  electron 
beams,  because  it  is  possible  to  carry  out  such  electric  discharge  uniformly 
at  3.5  atmospheres  of  pressure.  However,  the  electric  discharge  can  never  be 
said  to  be  practical.  This  is  because  the  operating  cost  is  high  and  a  larger 
device  must  be  used  for  the  electric  discharge.  Accordingly,  the  authors  have 
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begun  conducting  experiments  for  the  purpose  of  obtaining  greater  efficiency, 
higher  output,  and  high  repetition  by  using  a  more  practical  ultraviolet  [UV] 
preionization  system,  with  an  eye  toward  industrial  and  military  applications. 

This  report,  as  an  interim  report,  describes  output  characteristics  in  the 
case  where  the  total  pressure  in  a  laser  tube,  the  composition  ratio  expressed 
in  terms  of  a  partial  pressure  ratio  of  the  laser  gas,  the  reflectivity  of  the 
output  mirror,  and  the  charging  voltage  are  changed,  respectively,  because 
we  have  obtained  uniform  electric  discharge  in  the  atmosphere  and  have 
succeeded  in  obtaining  oscillation. 

2 .  Research  Method 

Figure  1  shows  schematic  drawing  of  an  experimental  device.  The  laser  tube  is 
made  of  transparent  acrylic  resin  and  is  sealed  off.  It  has  a  preionization 
pin  (cap  interval;  2.8  mm)  used  for  arc  discharge  and  a  main  electrode  made 
from  two  aluminum  bars  420  mm  long  and  20  mm  wide.  The  electrode  interval, 
discharge  length,  and  discharge  width  are  20,  390,  and  5  mm,  respectively. 


SG:  High  voltage  switch  =  SMfi,  Rj  =10kfl 

Cl  =  60  nF:  Main  condenser 

C2  =  1,700  pF  X  32:  Peaking  condenser 

This  device  has  a  capacity  transfer-type  excitation  circuit.  The  operation  of 
this  circuit  is  as  follows.  First,  the  main  condenser  (Ci — 60  nanofarads  (nF)) 
is  charged  from  a  high  voltage  power  source  through  a  current  limiting 
resistance  R^.  R2  is  a  charging  resistance.  Next,  when  the  spark  gap  (SG)  is 
switched  in  accordance  with  trigger  signals  given  from  a  trigger  circuit,  the 
electric  charge  from  will  charge  the  peaking  condenser  (C2)  of  1,700 
picofarads  (pF)  x  32  (8x2  rows  x  2  stages)  =  54  nF  through  the  preionization 
pin.  At  this  time,  UV  rays  will  be  emitted  from  the  preionization  pin,  charged 
particles  are  produced  between  the  main  electrodes,  and  preionization  occurs. 
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Simultaneously,  the  main  electrodes  begin  to  discharge,  energies  are  supplied 
from  Cl  and  C2,  and  laser  active  media  are  obtained. 

The  resonator  is  stable.  A  gold  mirror  made  of  Be-Cu,  with  a  curvature  of 
10  m,  was  used  as  a  totally  reflective  mirror,  while  a  plane  mirror  made  of 
fused  quartz,  with  a  thickness  of  10  mm,  was  used  as  the  output  mirror.  The 
diameter  for  both  mirrors  was  50  mm. 

The  GENTEC  ED-IOOA  was  used  as  a  highly  sensitive-joule  meter.  Laser  beams 
20  X  5  mm  were  focused  with  a  CaF2  lens,  and  measurements  were  carried  out, 
because  the  laser  output  is  comparatively  small.  The  total  pressure  in  the 
laser  tube,  the  composition  ratio  expressed  in  terms  of  the  partial  pressure 
ratio  of  the  laser  gas,  the  reflectivity  of  the  output  mirror,  and  the 
charging  voltage  of  Cj  were  all  changed  in  this  experiment.  The  laser  output 
change  in  the  respective  cases  was  measured  in  each  instance. 

3 .  Results 

The  experimental  results  are  shown  in  Figures  2  to  5.  Figure  2  shows  the  laser 
output  for  the  case  where  the  charging  voltage  is  changed,  using  total 
pressure  in  the  laser  tube  as  a  parameter.  Figure  3  shows  the  efficiency  for 
the  case.  Both  laser  output  and  efficiency  increase  as  total  pressure 
increases.  The  laser  output  has  the  maximum  value  for  the  charging  voltage, 
except  for  a  total  pressure  of  200  Torr,  and  efficiency  has  its  maximum  value 
in  the  same  way.  When  the  total  pressure  is  200  Torr,  even  if  the  charging 
voltage  increases,  the  laser  output  will  not  increase  but  in  fact  decreases, 
and  efficiency  also  decreases.  However,  it  is  possible  to  reduce  input 
energies  and  to  obtain  an  efficiency  that  exceeds  a  maximum  value  of  400  Torr, 
because  the  initial  discharge  voltage  is  low.  This  result  demonstrates  that  if 
input  energies  are  identical,  in  order  to  obtain  a  high  level  of  output  and  a 
high  degree  of  efficiency,  it  is  better  to  enhance  the  total  pressure. 
Further,  it  is  possible  to  obtain  the  charging  voltage  that  maximizes  laser 
output  by  means  of  this  total  pressure. 

Figure  4  shows  the  effects  of  changing  the  composition  ratio.  The  experiment 
was  conducted  by  using  two  kinds  of  partial  pressure  Ar:Xe  ratios — 100:1  and 
40:1.  This  experiment  demonstrated  that  when  the  amount  of  Xe  was  small,  the 
laser  output  was  large.  However,  it  is  believed  that  even  if  the  amount  of  Xe 
must  be  small,  an  optimum  value  exists,  because  when  there  is  no  Xe  there  will 
be  no  laser  output. 

Figure  5  shows  the  result  of  changing  the  reflectivity  of  the  output  mirror. 
In  this  case  the  charging  voltage  was  used  as  a  parameter.  It  was  learned 
that  an  optimum  value  exists  between  20—50  percent  in  both  cases. 

Figure  6  shows  an  example  of  oscillation  and  wavelength.  The  pulse  width  was 
about  1.2  US.  The  excimer  laser  has  a  pulse  width  of  several  tens  of  ns,^° 
while  this  laser  output  has  a  comparatively  long  oscillation. 
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OUTPUT  ENERGY  (mJ) 


VOLTAGE  (kV)  VOLTAGE  (kV) 

Figure  2.  Laser  Output  According  to  Figure  3.  Efficiency  According  to 
Charging  Voltage  When  Charging  Voltage  When 

Total  Pressure  Is  Regarded  Total  Pressure  Is  Regarded 

as  a  Parameter  as  a  Parameter 

Gas  composition  ratio  Ar:Xe  =  100:1 

Reflectivity  of  output  mirror  R  =  50  percent  (plane  mirror) 
Curvature  of  total  reflection  mirror  =  10  m  (gold  mirror) 


VOLTAGE  (kV) 

Figure  4.  Laser  Output  According  to  Charging  Voltage  When  Gas  Composition 
Ratio  Is  Regarded  as  a  Parameter 
Gas  composition  ratio  ArzXe 
Total  pressure  (TP)  =  760  Torr 

Reflectivity  of  output  mirror  R  =  50  percent  (plane  mirror) 
Curvature  of  total  reflection  mirror  =  10  m  (gold  mirror) 
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REFLECTIVITY  {%) 

Figure  5.  Laser  Output  According  to  Reflectivity  of  Output  Mirror  When 
Charging  Voltage  Is  Regarded  as  a  Parameter 
Gas  composition  ratio  Ar:Xe  =  100:1 
Total  pressure  (TP)  =  760  Torr 

Curvature  of  total  reflection  mirror  =  10  m  (gold  mirror) 


Figure  6 .  Example  of  Laser  Output  Waveform 

Gas  composition  ratio  Ar:Xe  =  139:1 
Total  pressure  TP  ^  700  Torr 

Reflectivity  of  output  mirror  R  =  50  percent  (plane  mirror) 
Curvature  of  total  reflection  mirror  =  10  m  (gold  mirror) 
Charging  voltage  V  =  30  kV 

4.  Discussion  and  Conclusion 

The  results  of  these  experiments  show  that  it  is  necessary  to  increase 
atmospheric  pressure  in  order  to  obtain  a  high  level  of  output  and  a  high 
degree  of  efficiency.  Also,  the  results  of  the  composition  ratio  experiment 
show  that  it  is  desirable  to  increase  the  amount  of  Ar  dramatically  rather 
than  that  of  Xe.  This  means  that  the  collision  pf  Ar  atoms  is  directly  related 
to  the  formation  of  the  laser  upper  order  and  that  it  is  difficult  to  form 
such  a  laser  upper  order  when  Xe  atoms  collide  with  each  other . ®  That  is ,  it 
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is  desirable  to  increase  the  niamber  of  Xe  atoms  in  a  state  where  only  Ar  atoms 
exist  around  Xe  atoms.  In  order  to  form  this  state,  it  is  necessary  to 
increase  not  only  the  pressure,  but  also  the  composition  ratio  in  conjunction 
with  the  increased  pressure.  However,  when  the  pressure  in  the  laser  tube  is 
iricreased,  main  discharge  will  not  readily  occur,  while  discharge  between  the 
cathode  electrode  and  the  pin  electrode  for  preionization  will  be  facilitated. 
For  this  reason,  it  is  necessary  to  study  the  conditions  under  which  this 
discharge  occurs  and  the  shape  of  the  electrodes  during  preionization,  and  to 
improve  the  laser  tube  to  prevent  this  discharge  from  readily  occurring. 

Also,  it  was  learned  from  the  experiment  that  the  accumulated  energies  (input 
energies)  are  comparatively  large,®-®  because  the  capacity  of  the  main 
condenser  is  60  nF.  These  energies  are  consumed  to  accelerate  ions  and 
electrons  ionized  in  the  preionization  system  and  to  further  ionize  many  Xe 
atoms  and  Ar  atoms.  It  is  considered,  however,  that  surplus  energies  formed 
after  these  Xe  and  Ar  atoms  are  ionized  may  cause  reionization  of  the  Xe  atoms 
excited  in  the  laser  upper  order.  It  is  also  believed  that,  because  of  this 
fact,  the  number  of  such  Xe  atoms  will  decrease  and  that  efficient  pumping 
will  become  impossible.  Accordingly,  in  order  to  form  the  laser  upper  order 
efficiently,  it  is  probably  necessary  to  form  optimum  kinetic  energies  for 
electrons . 

In  the  future,  we  plan  to  optimize  the  excitation  circuit,  improve  the 
preionization  system  with  a  view  to  obtaining  uniform  discharge  at  high 
atmospheric  pressure,  and  increase  the  laser's  output  and  efficiency  following 
the  recommendations  discussed  above. 
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Addition" ] 

[ Text ]  Summary 

The  burning  rate  of  AP  (ammonium  perchlorate)  composite  propellants  is 
decreased  significantly  by  the  addition  of  SrCOg  without  increasing  the 
pressure  exponent  of  the  propellant  burning  rate.  Results  of  microthermocouple 
measurements  show  that  surface  temperature,  Ts,  is  increased  from  700-970  K  by 
the  addition  of  SrCOg,  although  the  heat  release  from  the  gas  phase  to  the 
burning  surface,  (dT/dt)s+,  remained  unchanged.  DTA  (differential  thermal 
analysis)  and  TG  (thermo gravimetry)  experiments  were  also  carried  out  under  Ar 
flow  conditions.  The  AP  decomposition  temperature,  Tds ,  is  increased  from 
590-620  K  by  the  addition  of  SrC03. 

1.  Introduction 

The  rocket  motor  is  a  propulsion  engine  which  passes  gas  generated  by  the 
combustion  of  propellants  through  nozzles,  which  in  turn  discharge  the  gas  to 
the  outside  at  supersonic  speeds.  Thrust  is  produced  by  the  reaction 
generated  by  the  discharge.  A  particular  characteristic  of  the  solid  rocket 
motor  is  that  once  the  propellants  have  been  ignited,  it  cannot  be  controlled 
until  the  combustion  is  completed.  Therefore  the  gas  generating  speed  (i.e., 
burn  rate)  of  the  propellants  is  an  important  parameter  to  be  taken  into 
account  when  designing  a  rocket  motor. 

The  burn  rate  of  propellants  depends  on  their  physicochemical  properties,  and 
these  propellants  possess  limited  burning  rate  spectrum  characteristics. 
Therefore,  in  order  to  enhance  the  performance  of  rocket  motors,  it  is 
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necessary  to  increase  the  burn  area  of  propellants,  and  in  order  to  do  so,  it 
is  necessary  to  carry  out  research  on  increasing  and  decreasing  the  burn  rate 
of  propellants . ^ 

There  are  a  number  of  ways  to  increase  burn  rate.  These  include  increasing  the 
energy  contained  in  the  propellants ,  raising  the  thermal  conductivity  of  the 
burning  surface ,  and  chemically  promoting  the  burning  reaction  by  the  addition 
of  catalysts .  There  has  been  considerable  research  on  increasing  burn  rate, 
and  the  results  obtained  from  this  research  work  have  been  put  to  practical 
use,  taking  advantage  of  the  fact  that  this  increase  can  be  obtained  by 
raising  the  specific  impulse.^"® 

By  contrast,  little  research  on  decreasing  burn  rate  has  been  carried  out, 
because  it  is  very  difficult  to  lower  burn  rate  while  minimizing  decreases  in 
the  energies  contained  in  the  propellants . 

To  address  these  problems,  Glaskova,  et  al.,®  have  advocated  a  methodology  for 
lowering  the  burn  rate  by  controlling  the  chemical  reaction.  This  methodology 
indicates  the  following  methods  for  delaying  this  chemical  reaction  by 
altering  the  decomposition  of  ammonium  perchlorate  [AP] ,  which  is  an  oxidizer 
used  in  propellants . 

1)  Addition  of  a  substance  that  shifts  the  equilibrium  condition  of  the 
AP  decomposition  reaction  to  the  lift  side. 

2)  Addition  of  a  substance  that  combines  with  perchloric  acid  or  its 
decomposition  product  to  form  a  product  that  is  less  reactive  than 
perchloric  acid  or  AP. 

3)  Addition  of  a  substance  that  inhibits  the  decomposition  of 
perchloric  acid. 

The  use  of  these  substances  will  reduce  the  burn  rate  of  propellants,  because 
they  inhibit  the  oxidizing  reaction  of  ammonium. 

In  recent  years ,  there  has  been  a  tendency  for  the  operating  time  of  rocket 
motors  to  become  longer  (in  units  of  seconds).  For  this  reason,  it  is 
necessary  to  use  propellants  that  contain  greater  energies  than  conventional 
propellants  and  that  have  lower  burn  rates.  However,  there  has  been  little 
research  on  decreasing  burn  rate  through  the  use  of  so-called  negative 
catalysts . 

LiF  and  SrCOg  have  proven  to  be  effective  negative  catalysts  for  composite 
propellants  based  on  AP.  Hirata,  et  al.,^°  have  carried  out  research  on  this 
subject,  particularly  on  LiF,  and  have  identified  a  physicochemical  mechanism 
for  reducing  burn  rate.  However,  little  research  on  SrCOg  has  been  carried 
out,  and  a  similar  mechanism  has  not  been  identified.  In  our  research  work, 
we  have  identified  a  number  of  the  characteristics  of  SrCOg  and  have  studied 
its  effectiveness  in  reducing  burn  rate. 
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2 .  Sample  Propellant 


A  sample  propellant  was  made  on  a  trial  basis  to  investigate  the  influence  of 
strontium  carbonate.  This  work  was  undertaken  in  cooperation  with  the  Asahi 
Chemical  Industry  Co.  The  composition  of  this  propellant  is  shown  in  Table  1. 
Composition  N  is  for  reference  and  does  not  contain  any  strontium  carbonate. 
Composition  C  contains  2  percent  strontium  carbonate;  all  other  components  are 
the  same  as  in  composition  N. 


Table  1.  Composition  of  Sample  Propellant 


Propellant 

HTPB 

AP 

Diameter 

SrCOs 

400 

(/im) 

200 

10 

N 

12 

88 

34 

33 

33 

0 

C 

12 

88 

34 

33 

33 

2 

3 .  Experimental  Method 

3 . 1  Measurement  of  Burn  Rate 

The  sample  propellant  was  formed  into  a  long  strand  measuring  7  x  7  x  80  mm. 
It  was  coated  with  a  restrictor  made  of  silicon  rubber.  The  burn  rate  was 
measured  using  a  fuse  wire  method.  The  interval  of  the  fuse  wire  is  20  mm. 
The  propellant  was  pressurized  up  to  a  specified  pressure  by  injecting  N2  gas 
into  a  chimney- type  strand  combustor. 

3.2  Measurement  of  Burn  Rate  and  Temperature  Sensitivity 

The  experiment  was  conducted  by  setting  a  chimney-type  strand  combustor  in  a 
constant  temperature  bath  and  by  setting  the  strand  temperature  and  the  N2  gas 
temperature,  and  the  environmental  temperature,  at  243  and  343  K, 
respectively.  The  temperature  adjusting  time  was  set  at  72  hours  or  more! 
Also,  the  method  used  to  measure  the  burn  rate  is  the  same  as  that  described 
in  the  previous  item. 

3.3  Measurement  of  Combustion  Waves  by  Microthermocouple 

A  strand  testpiece  measuring  7  x  7  x  20  mm  was  divided  longitudinally  into 
halves  and  a  microthermocouple  with  a  diameter  of  12.5  /im  (Pt-Pt  10  percent 
Rh)  was  put  and  fixed  between  the  cut  surfaces  by  means  of  a  strand.  Poly¬ 
vinyl  alcohol  [PVA]  was  used  to  bond  the  strand.  The  strand  testpiece  was  left 
as  it  was  at  room  temperature  for  7  days  until  the  PVA  hardened  completely. 

It  is  possible  to  handle  homogeneous  propellants  like  double-base  propellants, 
because  their  flame  structures  are  homogeneous  and  can  be  regarded  as  one- 
dimensional  structures.  If  very  thin  thermocouples  are  used  in  such 
propellants,  the  measurement  results  obtained  from  time  areas  can  easily  be 
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converted  into  those  obtained  from  spatial  areas .  However ,  the  flame  structure 
of  composite  propellants  is  two-dimensional,  and  such  conversion  cannot  be 
done  readily.  This  is  because  these  composite  propellants  are  characterized 
by  a  diffusion  flame  generated  by  the  oxidizers  and  fuel  components.  Also,  the 
spatial  resolution  capability  is  insufficient,  because  both  the  particles  of 
the  oxidizers  and  the  thermocouple  junctions  have  limited  sizes.  A  qualitative 
assessment  was  made  on  the  basis  of  the  averaged  results  in  terms  of  time  and 
space  by  using  a  thermocouple  with  a  slightly  large  wire  diameter. 

3.4  Thermal  Analyzing  Experiment 

Differential  thermal  analysis  [DTA]  and  thermogravimetry  [TG]  were  carried  out 
simultaneously  as  thermal  analyzing  experiments  in  this  research  work.  The 
rate  at  which  temperature  rose  in  the  Ar  gas  atmosphere  (40  mi/min)  was 
20  K/min,  and  the  temperature  range  was  from  293-773  K.  The  A1  sealing  cell 
was  used  as  a  sample  cell. 

4.  Experimental  Results  and  Discussion 

4.1  Bum  Rate  Characteristics 

Figure  1  shows  temperature  sensitivity  characteristics  and  burn  rate 
characteristics  of  base  propellants  and  propellants  to  which  negative 
catalysts  were  added. 


Pressure,  p  (MPa) 

Figure  1.  Burn  Rate  Characteristics 

In  the  case  of  the  base  propellants ,  logarithms  of  the  burn  rate  were  changed 
almost  linearly  against  logarithms  of  pressure.  The  addition  of  SrC03  to  these 
base  propellants  will  lower  the  burn  rate. 
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Empirically,  the  burn  rate  of  propellants  is  represented  by  the  following 
equation: 

r  =  a  ■  p"  (1) 

Where , 

r :  Burn  rate 

a:  Constant 
p:  Pressure 

n:  Pressure  exponent  of  propellant  burn  rate 

Paying  attention  to  the  effectiveness  of  catalysts,  it  can  be  recognized  that 
the  effectiveness  of  a  decrease  in  burn  rate  is  remarkable  on  the  higher 
pressure  side.  That  is,  the  pressure  exponent  n  of  the  burn  rate  of  base 
propellants  is  lowered  from  0.77  to  0.50  at  a  temperature  of  343  K  and  from 
0.70  to  0.55  at  243  K,  respectively. 

LiF  has  proven  to  be  effective  as  a  negative  catalyst  affecting  the  burn  rate 
of  composite  propellants  based  on  AP.  It  is  known  that  the  addition  of  LiF  to 
such  propellants  is  unfavorable  in  practical  terms  because  adding  it 
dramatically  increases^®  the  pressure  exponent,  n. 

SrCOg  has  the  same  negative  catalyst  effectiveness  as  LiF.  SrCOj  is  believed 
to  be  more  superior  as  a  negative  catalyst  than  LiF  in  terms  of  decreasing  the 
pressure  exponent,  n. 

4.2  Temperature  Sensitivity  Characteristics 

Temperature  sensitivity,  Up  of  the  burn  rate  of  propellants  is  defined  by  the 
following  equations: 

CTp  ^  d  (In  r)/dTo  (2) 

Substituting  equation  (1)  for  equation  (2) , 

Op  =  1/a  •  (da/dTo)p  (3) 

Where , 

<Tp:  Temperature  sensitivity  of  burn  rate 

Tq:  Initial  temperature  of  the  propellant 

When  no  catalyst  is  added  to  a  propellant,  there  is  a  tendency  for  the 
temperature  sensitivity  to  rise  slightly  together  with  an  increase  in 
pressure,  but  when  negative  catalysts  are  added  to  propellants,  these 
propellants  will  lose  the  pressure  dependency  of  the  temperature  sensitivity, 
Op,  and,  as  shown  below,  the  temperature  sensitivity  will  become  almost 
constant,  as  shown  by  the  equation  ap  =  1.33  x  10“®  (1/K)  . 
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4.3  Study  of  Bum  Rate  Equation 

Analyzing  the  combustion  of  propellants  using  a  coordinate  system  based  on  the 
burning  surface  reveals  that  the  following  temperature  distribution 
(combustion  wave  structure)  is  formed. Figure  2(1)  is  an  example  of  a 
measurement  of  the  combustion  wave  structure  of  sample  propellants,  while 
Figure  2(b)  shows  a  typical  pattern. 


Figure  2(a).  Sample  Measurement  Figure  2(b).  Structural  Model  of 

of  Combustion  Waves  Combustion  Wave 

Using  a  Thermocouple 


The  initial  temperature,  Tq,  of  the  propellant  rises  sharply  in  the  vicinity 
of  the  burning  surface,  and  reaches  Tg  at  the  burning  surface.  Decomposition 
gases  generated  from  the  burning  surface  cause  an  exothermic  reaction  while 
flowing  downstream  at  a  certain  speed.  Therefore,  a  sharp  temperature  gradient 
is  formed  in  the  gas  phase,  and  the  temperature  reaches  the  final  flame 
temperature,  Tf .  When  part  of  the  heat  generated  from  the  gas  phase  is  fed  back 
to  the  burning  surface,  and  when  a  stable  combustion  wave  structure  is  formed, 
propellants  will  burn  at  a  steady  rate.  The  thermal  equilibrium  at  that  time 
is  represented  by  the  following  equation: 

Ap(dT/dx)g_  =  Ag(dT/dx)s+ 

+  Pp-r-Qs 

+  If  (4) 

Where , 

T :  Temperature 

x:  Distance 

r :  Burn  rate 

Qgi  Heat  of  reaction  at  the  burning  surface 

If:  Heat  feedback  from  the  luminous  flame  zone  to  the  burning 
surface  by  radiation 
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X:  Density 

p :  Thermal  conductivity 
p :  Propellant 

g:  Gas  phase 

S-;  Propellant  at  the  burning  surface 
s+:  Gas  phase  at  the  burning  surface 

The  left  side  of  equation  (4)  shows  the  thermal  movement  from  the  burning 
surface  to  the  solid  phase.  The  first  term  of  the  right  side  shows  the  thermal 
movement  from  the  gas  phase  to  the  burning  surface.  The  second  term  shows  the 
exotherm  or  the  endotherm  on  the  burning  surface,  and  the  third  term  shows  the 
thermal  movement  by  radiation  from  the  final  flame  to  the  burning  surface.  It 
is  generally  believed  that  the  third  term  of  the  right  side  can  be  ignored. 
The  combustion  wave  shown  on  the  left  side  is  two-dimensional,  and,  in  the 
case  of  steady— state  combustion,  the  following  equation  can  be  obtained  from 
the  law  of  the  conservation  of  energy: 

d(Xp- dT/dx)/dx  -  pp-r-dT/dx  =  0  (5) 

Integrating  equation  (5)  in  boundary  conditions  (x  =  0,  T  =  Tg,  and  x  =  -«i,  T 
=  Tg) ,  the  following  equation  can  be  obtained: 

Ap(dT/dx)s_  =  Ap-r-Cp(Ts  -  Tq)  (6) 

Where , 


c:  Specific  heat 

As  shown  in  the  following  equation,  the  burn  rate  can  be  represented  by 
combining  equations  (4)  and  (6): 

r  -  ^^idT/dx)^, 

PpCp(T^-To-QjCp)  ^  ^ 

Rewriting  this  equation,  the  following  equation  can  be  obtained: 


r  =  Qg  •  (8) 

Where , 

$  =  (dT/dx)g+  (9) 

^  =  (Tg-To-Qs/Cp)  (10) 

“s  “  ■^g/^p  ■  Pp)  (11) 


That  is,  the  causes  of  reduced  burn  rates  in  propellants  can  be  found  in 
increases  in  the  temperature  gradient,  (dT/dx)g+;  in  increases  in  surface 
temperature,  Tg  of  the  solid  phase;  or  in  decreases  in  the  exotherm,  Qg,  on  the 
solid  surface. 


32 


4.4  Change  of  Surface  Temperature 

The  surface  temperature,  Tg,  of  the  propellants  was  determined  on  the  basis  of 
the  results  of  the  strand  combustion  test  conducted  by  using  a  microthermo¬ 
couple.  Figure  3  shows  these  results. 


N 


o 

Figure  3 .  Combustion  Surface 


Temperature 


0.2  0.5  1.0  2.0 
Pressure,  p  (MPa) 

Figure  4.  Temperature  Change  Rate 
in  Gaseous  Phase 


Tg  of  base  propellants  does  not  depend  on  pressure  but  is  constant,  being  about 
700  K.  By  contrast,  it  was  ascertained  that  there  was  no  pressure  dependence 
in  propellants  to  which  negative  catalysts  were  added,  but  that  the  Tg  of  these 
propellants  was  increased  up  to  about  970  K.  It  is  believed  that  the  surface 
temperature  of  the  sample  propellants  rises  and  lowers  the  burn  rate. 


4.5  Change  of  (dT/dt)s+  by  Addition  of  Negative  Catalyst,  SrCOg 

As  shown  in  the  following  equation,  the  temperature  gradient  (dT/dx)g+  in  the 
gas  phase  can  be  represented  from  equation  (9) : 

dT/dx  =  dT/dt-dt/dx  =  r'^-dT/dt  (12) 

Therefore,  the  change  of  gas  phase  reactions  can  be  studied  by  using  the 
rising  temperature  rate,  (dT/dt)s+,  of  the  gas  phase  adjacent  to  the  burning 
surface  as  a  parameter.  (dT/dt)s+  can  be  determined  on  the  basis  of  the 
experimental  results  obtained  by  using  a  microthermocouple. 

Figure  4  shows  the  result  obtained  by  plotting  the  change  of  (dT/dt)g+  against 
pressure .  No  large  change  can  be  seen  in  spite  of  the  decrease  in  burn  rate . 
That  is,  it  is  believed  that  strontium  carbonate  does  not  affect  the  gas  phase 
reactions  of  propellants. 
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4.6  Exotherm,  Qg,  on  Burning  Surface 

Substituting  equation  (12)  for  equation  (7)  and  solving  for  Qg,  the  following 
equation  can  be  obtained: 

Qs  =  Cp-(Tg-To)  -  Ag/pp.(dT/dt)-(l/r)2  (13) 

Figure  5  shows  the  result  of  calculating  Qg  on  the  basis  of  equation  (13)  and 
the  measurement  results  that  have  been  obtained  up  to  now.  It  can  be  seen  that 
the  addition  of  SrC03  will  increase  about  40  percent  of  Qg  on  the  average . 
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Figure  5.  Heat  of  Reaction  of  Combustion  Surface 

4.7  Discussion  on  Negative  Catalyst  Effectiveness  for  Burn  Rate 

As  a  result  of  studying  the  combustion  wave  structure  adjacent  to  the  burning 
surface  on  the  basis  of  the  experimental  results  that  have  been  obtained  up  to 
now,  the  following  conclusions  can  be  drawn. 

Assuming  that  the  addition  of  catalysts  will  change  the  burn  rate  from  r^  to 
r^,  the  following  equations  can  be  obtained  by  using  equation  (8): 

=  “s.n  •  ^n/V>n  (14) 

=  «s,c  •  ^c/'Pc  (15) 

Where,  index  n  means  no  addition  of  catalysts,  and  index  c  means  addition  of 
catalysts.  In  the  same  way,  the  temperature  rising  rate  (dT/dt)g+  on  the 
burning  surface  can  be  given  by  the  following  equations  on  the  basis  of 
equation  (12): 


(dT/dt)g+_n  =  rn  ■ 

(dT/dx)g.^,n 

(16) 

(dT/dx)g+,c 

(17) 

Now,  assuming  that  the  percentage  of  change  in  the  burn  rate  obtained  by  the 
addition  of  catalysts  is  H,  the  following  equation  can  be  obtained  on  the 
basis  of  the  relations  in  equations  (8)  to  (11)  and  (14)  to  (17): 
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rc/^n 


N 


iT^.n-To-Os.n/Cr>)  (dT/ dt) 


(18) 


Let  the  physical  property  values  of  the  propellants  and  their  combustion  gas 
remain  unchanged  even  by  the  addition  of  catalysts  thereto.  According  to  the 
results  of  the  thermochemical  measurements  carried  out  in  this  research  work, 
it  can  be  understood  that  the  addition  of  SrCOg  will  have  little  effect  on  gas 
phase  reactions.  That  is,  the  following  equation  can  be  obtained,  because  the 
addition  of  catalysts  will  not  change  (dT/dt)s+,  which  shows  the  reaction  rate 
in  the  gas  phase: 


(dT/dt)s+,n  =  (dT/dt)s+.c 


(19) 


Therefore,  equation  (18)  can  be  represented  by  the  following  equation: 


= 


Ts.n-T, 


(20) 


The  heat  of  the  reaction  on  the  burning  surface  can  be  studied  by  using 
equation  (2)  on  the  basis  of  the  relationship  between  the  burn  rate  shown  in 
Figure  1  and  the  burning  surface  as  shown  in  Figure  4.  As  shown  in  Figure  5, 
the  addition  of  catalysts  will  increase  Qg  by  an  average  of  about  40  percent. 
Figure  6,  meanwhile,  shows  the  result  of  calculating  Qg  on  the  basis  of 
equation  (20). 


Pressure,  p  (MPa) 

Figure  6 .  Change  in  Burn  Rate 

There  is  a  tendency  for  an  increase  in  the  exotherm,  Qg,  on  the  burning  surface 
to  raise  the  burn  rate  as  shown  in  equation  (8).  However,  when  the  balance  of 
heat  adjacent  to  the  surface  of  the  propellants  is  considered  as  a  whole, 
(dT/dt)s  does  not  change,  but  Ts,c  becomes  sufficiently  large.  As  a  result,  E<1 
is  obtained  from  equation  (20),  and  the  burn  rate  is  lowered. 
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On  the  basis  of  the  above  discussion,  it  can  be  seen  that  an  increase  in  the 
burning  surface  temperatures  of  composite  propellants  based  on  AP  will  enable 
strontium  carbonate  (SrCOg)  to  reduce  the  burn  rate. 

4.8  Influence  of  Negative  Catalyst  on  Decomposition  Reaction  of  Oxidizer 

It  is  believed  that  the  addition  of  strontium  carbonate  will  affect  the 
decomposition  reaction  of  AP,  which  is  an  important  component  of 
propellants, 5'^''  Therefore,  these  mixtures  were  thermally  analyzed. 


Figure  7 .  Result  of  Thermochemical  Analysis 

Figure  7  shows  the  DTA  curves  and  TG  curves  of  strontium  carbonate,  AP,  and 
mixtures  of  these  elements.  Figure  7(a)  shows  the  result  of  a  thermal  analysis 
of  strontium  carbonate  as  a  simple  substance.  This  shows  that  no  change  is 
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observed  in  either  the  DTA  or  TG  curves  at  temperatures  ranging  up  to  730  K. 
However,  when  the  temperature  is  raised  to  1,280  K  by  using  an  open  cell  made 
of  aluminum  (AI2O3)  as  a  sample  cell,  an  endothermic  decomposition  reaction 
will  occur  at  temperatures  exceeding  1,000  K.  Calculating  the  weight  reducing 
rate  at  this  point,  the  following  equation  can  be  obtained: 

SrCOa  SrO  +  COaCg)  t  (21) 

This  result  corresponds  with  the  mass  of  CO2  in  the  case  where  such  a  reaction 
is  assumed.  Accordingly,  it  can  be  considered  that  the  above  reaction  occurs 
in  strontium  carbonate  as  a  simple  substance  at  temperatures  of  1,000  K  or 
more . 

Figure  7(b)  shows  the  result  of  a  thermal  analysis  of  AP  as  a  simple 
substance.  An  endothermic  peak  can  be  seen  in  the  vicinity  of  480  K  in  the 
DTA  curve  of  AP,  an  exothermic  reaction  begins  in  the  vicinity  of  590  K,  and 
a  large  exothermic  peak  can  be  seen  in  the  vicinity  of  740  K.  Looking  at  the 
TG  curve  in  this  case,  there  is  no  change  in  weight  at  the  endothermic  peak  at 
a  temperature  of  480  K.  The  endothermic  peak  at  this  temperature  is  due  to  the 
crystal  transformation  of  AP. 

In  comparing  the  DTA  and  the  TG  curves,  it  can  be  seen  that  both  weight 
reductions  start  simultaneously  at  a  temperature  of  590  K,  which  in  turn, 
starts  an  exothermic  reaction.  The  weight  reduction  at  the  exothermic  peak  is 
remarkable,  and  represents  99  percent  or  more  of  the  initial  loss  of  mass  of 
the  sample.  Considering  the  level  of  experimental  accuracy,  it  can  be 
considered  that  the  sample  decomposed  completely.  It  is  known  that  the 
decomposition  process  of  AP  is  a  multistage  reaction,  and,  finally,  all 
components  are  gasified  through  a  reaction  shown  in  the  following  equation. 
For  this  reason,  it  is  believed  that  the  results  shown  in  the  TG  curve  were 
obtained. 


NH4CIO4  -♦  NH3  +  HCIO4  (22) 

This  matter  is  the  result  of  thermally  analyzing  strontium  carbonate  and  AP  as 
simple  substances.  The  result  shown  in  Figure  7(c)  was  obtained  from  the 
thermal  analysis  curves  of  mixtures  made  by  adding  2  percent  (wt)  of  strontium 
carbonate  to  AP.  In  the  same  way  as  in  the  case  of  AP,  which  is  a  simple 
substance,  the  exothermic  peak  shows  a  crystal  transformation  point,  and  the 
temperature  in  this  case  shows  little  change. 

The  large  influence  of  the  addition  of  strontium  carbonate  can  be  seen  in  the 
temperature  range  where  an  exothermic  reaction  occurs.  A  slight  exothermic 
reaction  occurs  in  AP  when  used  as  a  simple  substance,  but  the  addition  of 
catalysts  will  cause  a  sharply  exothermic  reaction.  The  TG  curve  in  this  case 
shows  that  the  weight  reduction  begins  in  the  vicinity  of  520  K;  when  AP  was 
used  as  a  simple  substance,  this  onset  temperature  is  shifted  by  about  30  K  to 
the  high  temperature  side,  and  the  curve  inclination  is  large.  That  is,  the 
starting  time  and  the  weight  reducing  reaction  is  delayed  considerably,  and 
the  effectiveness  of  an  increase  in  this  reaction  rate  is  recognized. 
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Also,  the  exothermic  peak,  which  shows  the  maximum  exotherm,  occurs  at  about 
730  K.  Compared  with  AP  used  as  a  simple  substance,  this  temperature  is 
shifted  by  about  10  K  to  the  low  temperature  side.  Also,  the  exothermic 
research  is  particularly  sharp  in  this  case.  Seeing  the  Tg  curve  after  this 
exothermic  peak,  the  weight  reducing  rate  is  about  98  percent,  and  about  2 
Percent  remains  in  some  compound  form.  A  further  rise  in  temperatures  will 
gasify  this  substance. 


The  total  exotherm  caused  by  reactions  can  be  estimated  by  integrating  the  DTA 
curve.  Then,  the  respective  integrated  values  shown  in  Figures  7(b)  and  (c) 
were  compared  with  each  other,  but  little  change  was  seen  in  these  values. 

fact  indicates  that  the  addition  of  strontium  carbonate  will  not  change 
the  heat  insulating  flame  temperature  of  propellants. 

We  have  already  mentioned  that  the  weight  reduction  reaction  starting  point  is 
dramatically  shifted  to  the  high  temperature  side  in  the  TG  curve.  That  is, 
it  is  considered  that  strontium  carbonate  delays  the  start  of  the  exothermic 
decomposition  reaction  of  AP.  This  reaction  starting  point  temperature  was 
taken  as  Tjjg,  and  was  studied. 


1.5  1.6  1.7  1.8 

1000/T  (1/K) 

Figure  8.  Change  of  Activation  Energies 

Let  the  reaction  to  be  studied  be  an  Arrhenius-type  reaction.  Now,  the 
temperature  rising  rate,  6,  is  lower  than  the  burn  rate  of  pellets  in  the 
thermal  analyzing  test.  Figure  8  shows  the  result  of  studying  an  area  with  a 
higher  temperature  rising  rate  by  using  the  relation  {9  vs  1,000/T)  between 
the  temperature  rising  rate  {9)  and  actual  temperatures.  The  linear 
inclination  shown  in  Figure  8  indicates  activation  energies.  Therefore,  the 
activation  energy  of  AP  at  T^s  is  220  kJ/mole,  and  when  strontium  carbonate  is 
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added  to  AP,  the  figure  becomes  101  kJ/mole.  It  can  be  appreciated  that  the 
addition  of  strontium  carbonate  to  AP  will  lower  the  activation  energy  of  AP 
at  T(jg . 


When  the  straight  lines  shown  in  Figure  8  are  extended  up  to  a  temperature 
rising  rate  equivalent  to  the  burn  rate  of  propellants,  the  straight  lines 
shown  in  Figure  9  can  be  obtained.  Reading  the  temperatures  in  reverse  order 
from  the  respective  burn  rates  it  can  be  seen  that  the  temperature  in  the  case 
of  AP  as  a  simple  substance  is  about  1,100  K,  and  that  when  strontium 
carbonate  is  added  the  temperature  becomes  about  800  K.  This  result  accords 
almost  exactly  with  that  obtained  by  measuring  the  surface  temperature  by 
means  of  a  thermocouple.  That  is,  it  can  be  said  that  the  process  of  the 
reactions  generated  by  thermal  analyses  is  left  as  it  is,  and  the  reaction 
between  AP  and  strontium  carbonate  is  the  same  as  a  chemical  reaction  of 
burning  propellants . 

That  is,  the  addition  of  strontium  carbonate  to  AP  will  raise  the  decomposi¬ 
tion  temperature,  T^s,  of  AP.  As  a  result,  the  surface  temperature,  Tg,  of 
propellants  rises.  Subsequent  reaction  processes  are  accelerated,  because  the 
activation  energy  at  that  time  is  lowered.  As  a  result,  Qg  will  increase.  The 
reaction  of  the  gas  phase,  (dT/dt)g+,  is  measured  by  using  a  thermocouple,  but 
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there  is  no  change.  Therefore,  it  can  be  understood  that  strontium  carbonate 
does  not  affect  this  reaction.  Eventually,  the  change  in  the  decomposition 
process  of  AP  will  change  the  burn  rate  of  propellants. 

5.  Conclusion 

As  the  result  of  measuring  temperatures  by  means  of  a  thermocouple,  it  was 
ascertained  that  when  strontium  carbonate  was  added  to  composite  propellants, 
the  heat  of  the  reaction,  Qg,  and  the  surface  temperature,  Tg,  would  sharply 
increase  on  the  burning  surface,  and  the  temperature  rising  rate,  (dT/dt)g+, 
would  not  change  on  the  burning  surface.  The  reaction  between  AP  and  strontium 
carbonate  was  analyzed  and  studied  thermally.  As  a  result,  it  was  determined 
that  strontium  carbonate  raised  the  decomposition  temperature,  T^g,  of  AP.  It 
can  be  considered  that  this  reaction  system  is  the  same  as  that  in  the  case  of 
burning  propellants . 

That  is,  the  addition  of  strontium  carbonate  will  raise  the  decomposition 
temperature  of  AP.  The  gas  phase  of  the  combustion  wave  structure  has  no 
exothermic  change ,  and  the  change  of  heating  the  reaction  on  the  burning 
surface  is  comparatively  small .  With  regard  to  the  heat  balance  in  this  case , 
an  equilibrium  condition  is  attained  by  a  decrease  in  the  backward  rate  on  the 
burning  surface,  i.e.,  the  burn  rate  of  propellants. 
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Transonic  Wind  Tunnel  Test  of  Wing  Flutter  on  High  Maneuvering  Aircraft" ] 

[ Text ]  Summary 

It  is  of  concern,  with  regard  to  the  wing  flutter  of  aircraft  in  flight  at 
transonic  speeds,  that  there  is  a  phenomenon  called  "transonic  dip,"  which 
suddenly  reduces  the  flutter  speed  at  a  Mach  number. 


This  phenomenon  is  known  to  occur  primarily  on  the  wings  of  transport-class 
aircraft,  which  have  high  aspect  ratios  and  comparatively  thick  airfoils,  and 
it  appears  more  strongly  on  supercritical  airfoils  than  on  conventional  ones. 

A  transonic  wind  tunnel  test  was  conducted  to  investigate  what  characteristics 
of  the  flutter  occurred  on  other  types  of  wings,  those  of  fighter-class 
aircraft  which  have  small  aspect  ratios  and  thin  airfoils.  The  results  showed 
that  this  type  of  wing  was  also  affected  by  the  "transonic  dip"  with  the  same 
feature  for  the  airfoil,  and  that  the  external  store  under  the  wing  clearly 
reduced  the  flutter  speed  through  no  appearance  of  the  "transonic  dip"  could 
be  obviously  pointed  out. 

1.  Introduction 

There  has  been  extensive  research  on  the  various  kinds  of  wing  flutter 
experienced  by  aircraft,  and  the  flutter  generating  mechanism  has  been 
clarified.  Compressibility  affects  such  flutters  in  accordance  with  increases 
in  flying  speed.  That  is,  the  dynamic  flutter  pressure  (speed)  sharply 
decreases  in  the  vicinity  of  a  certain  Mach  number.  This  phenomenon  is  the 
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so-called  "transonic  dip,"  which  occurs  in  the  transonic  range.  It  is 
generally  said  that  when  a  supercritical  airfoil  (hereafter  "SC  airfoil")  is 
used  as  an  airfoil  on  the  main  wing  of  an  aircraft,  the  degree  of  this 
transonic  dip  will  be  greater  than  when  conventional  airfoils  are  used  on  the 
main  wings .  Research  concerning  transonic  dip  on  the  main  wings  of  many 
transport-class  aircraft  has  been  carried  out,  but  there  have  been  no 
comparable  studies  on  the  main  wings  of  fighter-class  aircraft. 

The  Third  Research  Center  has  tackled  this  problem.  It  has  made  a  model  for 
studying  the  flutter  of  main  wings  having  a  wing  planform  and  a  wing  section 
suitable  for  fighters.  It  has  conducted  a  transonic  wind  tunnel  test  of  the 
wing  flutter  using  this  model  and  has  confirmed  the  flutter  characteristics 
both  at  transonic  speeds  and  in  the  case  where  a  load  is  mounted  on  the  lower 
surface  of  a  wing. 


Angle  of  attack  (degree) 

Flutter  speed  (m/s) 

Equivalent  flutter  speed  (m/s) 

Flutter  dynamic  pressure  (kg/m^) 

Equivalent  flutter  dynamic  pressure  (kg/m^) 

Flexural  rigidity  (kgm^) 

Torsional  rigidity  (kgm^) 

3.  Model  as  a  Specimen 

Figure  1  shows  a  model  (hereafter  called  the  "main  wing  model")  used  to  study 
the  flutter  of  main  wings  having  a  wing  planform  and  a  wing  section  suitable 
for  high  maneuvering  aircraft. 

In  other  words ,  this  main  wing  model  has  the  same  wing  planform  as  that  of  the 
main  wings  of  a  wind  tunnel  model  used  for  research  (the  second  kind  of  A 
research  in  the  Research  Center)  on  the  shape  of  fighters  carried  out  by  the 
First  Research  Laboratory  of  the  First  Airframe  Department  (now  the  First 
Research  Laboratory  of  the  Aircraft  Department)  in  FY  1983.  The  main  wing 
model  is  a  half— mounting  model  reduced  by  about  8  percent  of  an  assumed  actual 
aircraft  with  a  wingspan  of  about  10.5  m,  an  aspect  ratio  of  about  3.14,  and 
a  leading  edge  sweepback  angle  of  about  32.5  degrees.  Table  1  shows  the 
degree  of  similarity  between  the  model  as  a  specimen  and  hypothetical  actual 
aircraft . 

As  shown  below,  such  main  wing  models  were  manufactured  every  fiscal  year. 

Five  kinds  of  main  wing  models  were  manufactured  in  FY  1983.  They  have 
NACA65A004.8  Ma  airfoils,  as  shown  in  Table  2  and  Figure  2,  and  have  different 
rigidities . 


2 .  Symbols 


Vp: 

VpEAS- 

qp: 

qPEAS^ 

El: 

GJ: 
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Table  2 .  Wind  Tunnel  Model  Wing  by  Type  and  Coordinates 


X 


0 


3 


3  2.5  0 


3  7.5  0 


4  0.0  0 


4  2.5  0 


4  5.00 


5  2.5  0 


5  5.0  0 


5  7.5  0 


6  0.0  0 


6  2.5  0 


6  5.0  0 


6  7.5  0 


7  0.0  0 


7  2.5  0 


7  5.  0  0 


7  7.5  0 


8  0.0  0 
8  2.5  0 


8  5.00 


8  7.5  0 


9  0.0  0 


9  2.5  0 


9  5.0  0 


9  7.5  0 


Standard  airfoil 


Y\  ‘f>C 


0  0 


SC  airfoil 


Yu  5SC 


0 


Y1 


-0.2  607 


0.4  3  6  8 


0.4  9  7  9 


-0.5  9  6 

2 

-0.6  3  8 

IIIIK 

.6  75  8 

0.8  0  1  2 


.8  5  30 


— 

0.8  9  9  6 

— 

0.9  4  2  0 

0.9  8  0  8 

— 

miiiM 

— 

1.28  0  8 

— 

1.3  8  4  6 

.0 

0 

1  2.35  6  8  1 

1  -2.3  5  68 

i  3.2  6  9  2  1 

-1.7  3  8 

3 


6 


3.2  0  0  0 


3.1  5  38 


3.0  9  5  2 


3.0  2  6  9 


-1.7  30  8 


0.9  2  7  1 


•0.9  27  1 


1.0  5  00 


0.8  6  1  5 


0.6  6  9  2 


0.4  7  3  1 


23 


-1.4  80  5 


-1.4  0  77 


-1.32  31 


-1.1  192 


-1.0  0  77 


-0.8  92  3 
-  0.7  7  6  9 


5 


-0.5  4  62 


-0.4  3 


-0.32  31 


-0.10  00 


0.2  6  9  2 

-0.1  0  7  7 

0.0  5  7  7 

-0.196  2 

-0.1  65  4 

-0.3  6  15 

0.33  0  0 

45 


sc  airfoil 


I  I 

60 - 7  0' 


10(MtC 


Symmetrical 

airfoil  (  65A004.8Md) 


— 1 0 - 20 - 30 - 4  0 - 5  0 - 6  0- - 7  0 - 80"  9o'  10096C 

Figure  2.  Airfoil  of  Model  as  a  Specimen 

Five  kinds  of  main  wing  models  were  manufactured  in  FY  1984.  They  have  SC 
airfoils,  as  shown  in  Table  1  and  Figure  2.  Their  wing  thickness  is  almost 
the  same  as  that  of  the  main  wing  models  manufactured  in  the  previous  fiscal 
year,  and  have  different  rigidities. 

Three  kinds  of  main  wing  models  were  manufactured  in  FY  1985.  They  have 
NACA65A004.8  Md  airfoils,  and  have  different  rigidities.  In  addition,  a 
similar  scale  model  of  a  slip  tank  with  a  capacity  of  620  gallons  was  made  as 
a  load  mounted  on  the  lower  surface  of  the  actual  aircraft  by  using  a  pylon. 

Next,  an  aluminum  core  was  inserted  in  the  center  of  the  main  wing  model.  The 
mass  and  rigidity  distributions  of  the  wings  were  represented  by  this  aluminum 
core,  which  was  wrapped  and  formed  by  a  foamed  material,  and  eight  notches 
were  made  in  the  direction  of  the  wing  chord  line  of  the  wing  span  to  prevent 
the  foamed  material  from  becoming  rigid. 

Also,  as  previously  mentioned,  the  purpose  of  this  test  was  to  confirm  the 
influence  of  the  airfoil  or  the  load  mounted  on  the  lower  surface  of  the  main 
wing  or  transonic  flutter  characteristics.  Therefore,  the  wash-out  angle  and 
dihedral  angle  of  the  main  wing  model  were  zero  degrees.  This  model  also  has 
an  identical  airfoil  from  the  wing  root  and  wing  tip  to  avoid  other 
aerodynamic  influences . 

In  addition.  Figure  3  shows  a  model  of  a  slip  tank  as  a  load  mounted  on  the 
lower  surface  of  the  main  wing.  The  aluminum  core  was  wrapped  and  formed  with 
wooden  materials,  and  the  shape,  distribution  of  mass,  and  moment  of  inertia 
were  simulated  to  resemble  those  of  an  actual  aircraft. 
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(2 


Main  wing 

@  Slip  lank  (wooden  material:  a  heavy  weight  is  placed 
— .  inside) 

Pylon  (wooden  material)  (Unit:  m) 


Figure  3.  State  of  Slip  Tank  Mounted  on  Aircraft 
4.  Wind  Tunnel  Test 

The  wind  tunnel  test  was  conducted  in  the  high-speed  wind  tunnel  facility  of 
Fuji  Heavy  Industries'  [FHI]  Utsunomiya  Manufacturing  Division. 
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4 . 1  Measuring  Method 


A  strain  gauge  was  attached  to  the  wind  tunnel  wall  where  the  aluminum  core  of 
the  main  wing  model  was  attached.  The  output  was  measured  with  the  measuring 
system  shown  in  Figure  4.  The  power  spectra  of  vibration  were  analyzed,  and 
the  vibration  frequency  strain  and  test  states  were  investigated.  This 
experiment  made  it  possible  to  identify  the  generation  of  flutter. 


Model  as  a  specimen 


Figure  4.  Block  Diagram  of  Instrumentation  System 
4.2  Flutter  Test 

The  main  wing  model  shown  in  Figure  3  was  installed  on  a  dummy  fuselage  of  the 
wind  tunnel  wall,  and  the  angle  of  attack  was  set  to  a  =  0°.  The  flutter  test 
was  started  at  a  Mach  number  lower  than  that  at  which  the  generation  of 
flutter  was  predicted  by  analytical  calculation,  and  a  Mach  number  at  which 
such  generation  began  was  found  by  gradually  increasing  the  Mach  number. 

Further  increases  in  the  Mach  number  beyond  the  Mach  number  at  which  flutter 
appears  will  not  generate  flutter,  depending  on  the  wing  shape.  The  flutter 
test  was  conducted  at  around  Mach  1.0  to  confirm  the  existence  of  Mach  numbers 
for  which  the  generation  of  flutter  disappears  by  Mach  1.0.  When  no  flutter 
was  generated,  the  Mach  numbers  were  lowered  gradually  to  find  a  Mach  number 
at  which  flutter  appeared. 
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Also,  the  main  wing  model  was  excited  by  using  natural  disturbances  in  the  air 
flow  in  the  wind  tunnel,  and  when  the  slip  tank  was  mounted  on  the  lower 
surface  of  the  main  wing,  the  wing  tip  was  drawn  with  enamel  wires. 

Figure  5  [not  reproduced]  is  a  photograph  of  the  wind  tunnel  test. 

5.  Results  and  Study 

The  results  of  the  wind  tunnel  test  are  shown  in  Figure  6. 


Mach  number 


Figure  6.  Relationship  Between  Mach  Number  and  Flutter  Velocity 

Based  on  the  test  using  the  model  shown  in  Figure  5  [not  reproduced] ,  it  was 
confirmed  that  transonic  dip  (hereafter  called  "dip")  occurred  even  in  main 
wings  equipped  with  airfoils  having  a  small  aspect  ratio  and  a  small  wing 
thickness  (about  5  percent  C) . 

It  was  demonstrated  that  a  dip  of  about  5  percent  existed  in  main  wings 
equipped  with  conventional  airfoils  in  the  vicinity  of  M  =  0.87,  and  that  a 
dip  of  about  10  percent  existed  in  those  having  SC  airfoils.  However,  both 
Mach  numbers  that  show  the  existence  of  such  dips  seem  to  be  almost  identical. 
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Also,  as  generally  mentioned,  the  dip  in  the  case  of  the  main  wings  equipped 
with  SC  airfoils  is  larger  than  that  in  the  case  of  those  having  conventional 
airfoils.  It  should  be  noted,  however,  that  some  of  the  differences  in  the 
amount  of  dip  are  affected  by  the  presence  or  absence  of  cambers .  This  is 
because,  in  the  above  case,  the  conventional  airfoil  is  symmetrical,  while  the 
SC  airfoil  has  a  cambered  structure  as  well  as  SC  pressure  distribution. 

In  addition,  the  amount  of  dip  existing  in  the  main  wings  of  transport— class 
aircraft  is  said  to  be  20-30  percent,  while  for  this  model  the  value  was  as 
small  as  10-15  percent.  We  believe  that  this  is  because  the  aspect  ratio  and 
wing  thickness  of  the  main  wing  of  the  model  are  half  those  of  the  main  wing 
of  transport— class  aircraft. 

In  general,  there  is  a  tendency  for  the  flutter  speed  to  lower  when  a  slip 
tank  is  mounted  on  the  lower  surface  of  the  main  wing  as  opposed  to  the  case 
when  no  slip  tank  is  mounted  (hereafter  called  "clean  state").  This  tendency 
was  clarified  by  this  test,  but  the  amount  of  the  dip  was  not  clarified.  This 
is  because  the  tendency  was  confirmed  only  at  three  points  (mark  •  shown  in 
Figure  5  [not  reproduced]). 

In  addition,  a  comparison  with  the  estimated  value  was  made  using  nonsteady 
aerodynamic  force  calculated  by  the  doublet  lattice  method  [DLM] .  The 
estimated  flutter  value  (hereafter  called  "Value  Estimated  by  Modified  Strip 
Analysis  [MSA]")  was  calculated  by  PASPA.  A  comparison  between  the  test  value 
and  this  estimated  flutter  value  was  made  by  using  the  DLM  estimated  flutter 
value  calculated  by  the  NASTRAN  and  by  using  the  nonsteady  aerodynamic  force 
calculated  by  MSA  (Figure  6) . 

The  value  estimated  by  the  DLM  accords  well  with  the  tested  value  in  the  clean 
state  except  for  the  vicinity  of  the  dip.  By  contrast,  the  value  estimated  by 
using  MSA  is  considerably  lower  than  the  test  value. 

There  was,  however,  a  large  difference  between  the  value  estimated  by  using 
the  DLM  and  that  by  using  the  MSA  when  a  slip  tank  is  mounted  on  the  lower 
surface  of  the  main  wing.  That  is,  in  the  case  where  the  value  was  estimated 
by  the  DLM,  the  flutter  speed  was  higher  than  that  in  the  clean  state,  but  in 
the  case  where  the  value  was  estimated  by  the  MSA,  the  flutter  speed  was  lower 
than  that  in  the  clean  state.  The  value  estimated  by  the  DLM  was  higher  than 
that  for  the  clean  state.  As  a  result,  it  was  possible  to  make  the  following 
determination . 

Assuming  that  the  rigidity  is  constant,  the  flutter  speed  is  affected  greatly 
by  the  mass  and  moment  of  inertia  of  loads,  and  by  the  relative  position 
between  the  elastic  axis,  loading  axis,  and  center  of  gravity,  etc. 

Figures  7  to  9  show  the  results  of  parametrical  estimates  for  each  case  by 
both  DLM  and  MSA. 
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Figure  7.  Influence  of  Mass  Figure  8.  Influence  of  Moment  of 

Inertia 

(1)  When  the  additional  mass  is  increased  gradually  from  0,  both  DIM  and  MSA 
predict  that  the  flutter  speed  will  decrease,  subsequently,  the  flutter  speed 
will  increase  at  a  certain  point,  and  this  is  where  the  difference  between  the 
values  estimated  by  DIM  and  MSA  becomes  large. 

(2)  When  the  value  of  the  moment  of  inertia  is  0  at  a  constant  mass,  the 
estimated  values  for  DLM  and  MSA  will  be  close  to  each  other,  but  there  is  a 
tendency  for  the  DLM  estimated  value  to  sharply  increase  at  a  certain  point. 

(3)  When  the  center  of  gravity  of  loads  is  at  the  back  of  the  elastic  axis  of 
the  main  wing,  there  is  little  difference  between  the  estimated  values  for  DLM 
and  MSA,  but  when  it  is  moved  further  toward  the  front  of  the  main  wing,  there 
is  a  tendency  for  the  difference  in  estimated  values  to  become  large. 

An  examination  of  items  (1),  (2),  and  (3)  indicates  that  the  model  as  a 
specimen  provides  data  on  state  that  causes  the  difference  between  values 
estimated  by  using  both  DLM  and  MSA  becomes  unexpectedly  large.  It  is 
believed  that  this  difference  may  be  the  result  of  the  different  means  by 
which  DLM  and  MSA  calculate  nonsteady  aerodynamic  force,  but  no  clear  cause 
has  yet  been  determined. 
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<—  Front  Position  of  center  of  Rear 
gravity  of  slip  tank 
(m) 

Figure  9.  Influence  of  Position  of  Center  of  Gravity 
6.  Discussion 

In  this  wind  tunnel  test,  as  shown  in  item  5.1,  the  generation  of  flutter  was 
judged  by  analyzing  vibration  spectra  and  monitoring  the  state  and  the  point 
at  which  the  reciprocal  number  of  strain  per  hertz  of  vibration  raised  to  the 
second  power  sharply  approaches  0. 

It  was  confirmed  that  there  was  a  tendency  for  this  reciprocal  number  to 
approach  sharply  0  in  the  case  of  the  clean  state.  It  was  also  possible 
visually  to  identify  a  sign  whereby  the  model  as  a  specimen  causes  flutter. 
However,  when  a  slip  tank  was  mounted,  the  reciprocal  number  did  not  approach 
0.  Rather,  the  model  as  a  specimen  gradually  entered  a  state  of  permanent 
vibration,  and  there  was  no  tendency  toward  divergence.  Accordingly,  it  was 
difficult  to  recognize  the  generation  of  flutter.  It  is  believed  that  this  is 
because  the  wind  tunnel  blockage  ratio  is  about  1  percent  for  the  clean  state 
and  about  2  percent  when  the  slip  tank  was  mounted,  thus  air  flow  energies  in 
the  wind  tunnel  were  not  able  to  cause  divergent  vibrations. 
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Next,  the  model  has  a  simple  structure  and  is  economical,  but  has  a  fixed 
rigidity.  For  more  thorough  testing,  it  would  be  necessary  to  manufacture 
several  wings  with  different  rigidities.  Therefore,  the  goal  of  this  wind 
tunnel  testing  method  should  be  to  simulate  several  different  wings  completely 
and  aerodynamically.  Accordingly,  the  wind  tunnel  testing  method  would  be  able 
to  obtain  a  rough  congruence  between  the  dip  produced  by  flutter  generated 
from  actually  designed  wings  and  wings  used  as  the  objects  of  research  work. 

The  results  of  the  wind  tunnel  test  demonstrate  that  the  value  estimated  by 
using  DLM  accords  well  with  the  test  value,  except  for  the  portion  of  the  dip 
in  the  clean  state,  but  that  is  the  opposite  of  the  test  value  when  a  slip 
tank  is  mounted.  The  same  tendency  applies  to  the  values  estimated  by  using 
MSA,  but  it  is  necessary  to  pay  attention  to  the  fact  that  the  value  estimated 
by  using  the  DLM  is  higher  than  the  test  value  in  portions  which  seem  to  be 
dips  in  the  slip  tank  mounting  state.  Therefore,  it  would  be  permissible  to 
adopt  the  value  estimated  by  using  MSA  in  the  clean  state  in  order  to  ensure 
the  safety  of  actual  aircraft  against  transonic  flutter. 

7.  Conclusion 

We  have  made  a  flutter  wind  tunnel  test  model  (half— mounting  model)  with  an 
airfoil  and  a  wing  platform  suitable  for  high  maneuvering  aircraft,  have 
conducted  a  transonic  flutter  wind  tunnel  test,  and  have  made  a  comparison 
between  flutter  values  estimated  by  using  DLM  and  MSA,  and  actual  test  values. 
DLM  is  used  to  calculate  nonsteady  aerodynamic  forces  in  the  NASTRAN  program, 
while  MSA  is  used  to  calculate  nonsteady  aerodynamic  forces  in  the  PASPA 
program.  The  following  conclusions  were  drawn  from  the  results  of  the 
comparison. 

(1)  The  existence  of  dips  was  confirmed.  As  previously  noted,  it  was 
demonstrated  that  the  magnitude  of  dips  generated  from  SC  airfoils  was  larger 
than  that  from  conventional  airfoils.  However,  the  difference  between 
conventional  and  SC  airfoils  includes  the  difference  between  symmetrical 
airfoils  and  camber  ones. 

The  Mach  number  for  dips  generated  by  the  test  wing  was  in  the  vicinity  of 
M  “  0.87,  and  the  magnitude  of  dips  in  the  case  of  conventional  airfoils  was 
about  5  percent,  while  that  in  the  case  of  SC  airfoils  was  about  10  percent. 

(2)  With  regard  to  the  comparison  between  the  test  value  and  the  estimated 
value,  the  value  estimated  by  using  DLM  accorded  well  with  the  test  value  in 
the  clean  state,  except  in  the  vicinity  of  the  dip.  By  contrast,  the  value 
obtained  by  using  MSA  was  considerably  lower  than  the  test  value. 

The  value  estimated  by  using  DUt  when  a  slip  tank  was  mounted  was  higher  than 
that  in  the  clean  state,  and  the  value  estimated  by  MSA  when  a  slip  tank  was 
mounted  was  lower  than  that  in  the  clean  state.  This  tendency  generally 
accorded  with  the  actual  test  value,  but  appeared  to  be  higher  than  the  test 
value  in  the  portion  that  seemed  to  be  a  dip . 
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Appendix 


In  order  to  conduct  wind  tunnel  tests  while  arbitrarily  changing  the  flutter 
dynamic  pressure  against  Mach  numbers,  it  is  necessary  to  vary  the  density  of 
the  wind  tunnel  and  to  widen  the  range  of  the  variable  density. 

Further,  when  wind  tunnel  tests  are  conducted  by  using  a  scale  model  (this  is 
called  a  "completely  similar  model")  made  by  completely  simulating  the 
dimensions,  rigidity,  and  moment  of  inertia  of  the  main  wing  of  the  assumed 
actual  aircraft,  in  general,  no  flutter  is  generated  from  the  main  wing  of  the 
scale  model.  However,  in  order  to  obtain  actual  transonic  flutter  characteris¬ 
tic  curves  from  this  completely  similar  model,  wind  tunnel  tests  are  conducted 
by  using  the  above-mentioned  variable-density  wind  tunnel,  or  a  number  of 
models  to  be  used  as  specimens  are  manufactured  to  cause  flutter  by  lowering 
the  rigidity  of  the  completely  similar  model,  a  Mach  number  which  generates 
flutter  is  tentatively  found,  and  characteristic  curves  are  determined  by 
converting  this  Mach  number  for  application  with  the  completely  similar  model. 
The  latter  method  was  adopted  in  this  test.  That  is,  assuming  that  the  flutter 
dynamic  pressure  of  the  completely  similar  model  is  qp,  that  of  the  low 
rigidity  model  is  qptesti  the  estimated  flutter  dynamic  pressure  of  the 
completely  similar  model  is  qp^,  the  estimated  flutter  dynamic  pressure  of  the 
low  rigidity  model  is  qpLc,  and  the  Mach  number  at  that  time  is  Mp,  the 
following  equation  can  be  obtained: 

QF=-^'QFtest 

'^FLC 

The  flutter  dynamic  pressure  of  the  completely  similar  model  against  the  Mach 
number  Mp  obtained  from  the  test  can  be  obtained  from  this  equation. 


Figure.  Relationship  Between  Rigidity,  Wind  Tunnel 
Dynamic  Pressure,  and  Mach  Number 
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The  flutter  dynamic  pressure  is  regarded  as  a  wind  tunnel  dynamic  pressure, 
because  the  wind  tunnel  dynamic  pressure  can  be  obtained  directly  from  the 
wind  tunnel  test.  This  manuscript  describes  the  flutter  speed  as  flutter 
dynamic  pressure,  because  it  seems  to  be  easy  for  readers  to  understand  the 
concept  of  flutter  dynamic  pressure. 
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[Text]  Sxxmmary 

This  paper  describes  fundamental  studies  and  test  results  on  high-rate 
dischargeable  lithium-thionyl  chloride  batteries  fabricated  in  1985-1986. 

1.  Fundamental  Studies 

Temperature  dependence  on  electric  conductivity  of  several  electrolyte 
solutions  was  studied,  so  that  the  optimum  electrolyte  was  selected  for  the 
battery.  And,  furthermore,  we  established  prospects  of  preventing  voltage 
delay  at  initial  discharge  by  coated  negative  lithium  electrodes. 

2.  Test  Results 

Cylindrical  batteries  and  disk— shaped  batteries  were  fabricated.  The 
cylindrical  battery  focuses  increased  power  capacity  and  has  spiral  bobbin¬ 
shaped  electrodes.  The  test  results  show  the  capacity  of  10.3  Ah  and  31.5  Wh 
at  high  current  density  of  10  mA/cm^.  The  disk-shaped  battery  focused  improved 
resistance  against  shock  and  vibration  for  a  projection-type  wave  jammer.  The 
test  results  show  no  problems  under  vibration  of  10-500  Hz  (10  G)  .  In 
addition,  we  established  that  the  batteries  can  discharge  at  a  high  current 
density  of  10  mA/cm^  at  a  low  temperature  of  -40°C.  The  assembled  batteries, 
four  connected  disk-shaped  batteries  in  series,  were  fabricated  and  tested  for 
confirmation  of  safety. 
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1.  Introduction 


The  realization  of  primary  batteries  with  high  capacity  and  highly  efficient 
electric  discharge  is  necessary  to  cope  with  specific  power  sources  that  will 
be  used  in  the  future  together  with  the  miniaturization  of  electronic 
equipment  and  increases  in  performances . 

Research  on  lithium-thionyl  chloride  batteries  has  been  underway  since  FY 
1984.  First,  the  possibility  of  highly  efficient  electric  discharge  was 
demonstrated^  by  using  coin-type  batteries.  Next,  output  density  was  increased, 
and  mechanical  strength  and  safety  were  established  as  objects  of  this 
research  work.  With  regard  to  the  increase  in  output  density,  basic  electrical 
characteristics  were  studied  in  FY  1985  by  using  a  cylindrical-type  (unit  1) 
battery  with  a  spiral  electrode  structure  shown  as  in  Figure  1  and  Photograph 
1  [not  reproduced] .  With  regard  to  mechanical  strength,  vibration  and  impact 
resistances  were  confirmed  in  FY  1986  by  using  a  battery  with  the  disk-type 
electrode  structure  shown  in  Figure  2  and  Photograph  2  [not  reproduced] . 
Research  on  the  following  topics  was  carried  out  during  the  same  fiscal  year: 


Laser  welding 


No 

Name  of  parts 

Remarks 

1 

Battery  storing  can 

SUS304 

2 

Battery  cap 

Fe-Ni  alloy 

3 

Positive  electrode 

Carbon  rack 

4 

Negative  electrode 

Lithium 

5 

Separator 

Nonwoven  glass  fabric 

6 

Upper  portion  insulating  sheet 

tt  tt  It 

7 

Lower  portion  insulating  sheet 

If  II  II 

8 

Negative  electrode  terminal 

9 

Positive  electrode  tab 

11  tt  It 

10 

Negative  electrode  tab 

tt  tt  tt 

11 

Ceramic 

Alumina 

12 

Sealer 

Epoxy-resin 

Figure  1.  Cylindrical  Lithium-Thionyl  Chloride  Battery 
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fl  battery 


Dimension:  m 
I  Weight:  about  i9G  g 


±2.0 


76.0 


B  battery 


+  I 


1 

13.5±1 

Dimension: 

m 

Weight: 
about  27g 

......  ,n 

fl' 
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Laser  welding 


No 

Name  of  parts 

Remarks 

Bi 

Battery  storing  can 

SUS304 

2 

Positive  electrode 

Carbon 

Negative  electrode 

Lithium 

01 

Separator 

Nonwoven  glass  fabric 

5 

Insulating  ceramic 

Alumina 

6 

Negative  electrode  terminal 

Nickel 

Figure  2.  Disk-Type  Lithium-Thionyl  Chloride  Battery 


1)  a  voltage  delay  phenomenon  generated  at  the  initial  period  of  electric 
discharge,  and  2)  methods  of  preventing  this  phenomenon.  A  demonstration  test 
of  the  safety  was  conducted  in  FY  1987  by  using  a  12-V  practical  battery  like 
that  shown  in  Figure  3  and  Photograph  3  [not  reproduced] .  This  battery  is  an 
aggregation  of  unit  cell  batteries.  This  report  will  hereunder  describe  the 
above  research  processes  and  results. 


2 .  Experiment  and  Results 


2.1  Clarification  of  Elemental  Technologies 

(1)  Values  for  resistance  against  temperature  were  measured  and  electrical 
conductivity  was  determined  by  calculating  the  concentration  of  electrolyte 
(solvent;  SOCI2  and  solute:  LiAlCl^)  as  a  parameter.  The  results  are  shown 
in  Figure  4. 
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Outline  and  dimension  of  B  battery 

Figure  3.  Aggregated  Lithium-Thionyl  Chloride  Battery 
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Figure  4.  Temperature  Dependency  of  Electric  Conductivity 
(2)  Prevention  of  generation  of  voltage  delay  phenomenon 

The  voltage  delay  phenomenon  means  that  a  threshold  voltage  delay  or  voltage 
drop  is  generated  during  the  initial  period  of  electric  discharge  by  an  LiCl 
coat  formed  on  the  Li  surface  of  the  negative  electrode.  This  phenomenon 
occurs  primarily  in  the  following  cases: 

i)  The  storage  period  is  long. 

ii)  The  storage  temperature  is  high. 

iii)  The  discharge  current  density  is  large 

iv)  The  ambient  temperature  is  low. 

The  following  can  be  cited  as  methods  of  preventing  this  phenomenon: 

a)  The  Li  surface  of  the  negative  electrode  is  coated  so  it  will  not 
form  an  LiCl  coat. 

b)  An  electrolyte  based  on  Li20  +  AICI3  is  used. 

As  a  result  of  studying  these  two  methods,  it  was  recognized  that  both  methods 
improve  the  effectiveness  of  the  batteries. 

(Coating  Method) 

Two  kinds  of  batteries  were  trial  manufactured  as  specimens  by  using  cyano¬ 
acrylate  polymer  [CA]  and  polyvinyl  chloride  [PVC]. 

Figure  5  shows  the  initial  voltage  characteristics  for  the  electric  discharge 
of  two  batteries  in  which  the  Li  pole  of  one  was  not  treated  while  the  other 
was  coated  with  CA.  Both  batteries  were  stored  at  a  temperature  of  45°C  for 
3  days  and  then  were  discharged  electrically  at  room  temperature  at  a  current 
density  of  10  mA/cm^.  The  voltage  of  the  battery  coated  with  Ca  rose  in  a 
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Figure  5.  Initial  Voltage  Characteristics  of  Electric  Discharge 


moment  of  electric  discharge,  but  that  in  which  Li  pole  was  not  treated  had  a 
voltage  delay  time  of  3 . 3  seconds  until  the  voltage  recovered  to  2  V. 


[Electrolyte  Method] 

It  is  necessary  to  study  this  method  carefully,  because  there  are  problems  in 
the  discharge  voltage  and  discharge  duration. 

2.2  Experiment  Using  Test  Batteries 

2.2.1  Electric  Discharge  of  Cylindrical-Type  Battery 

A  cylindrical-type  battery  whose  electrode  has  a  spiral  bobbin  structure  was 
adopted  to  make  it  possible  to  increase  packaging  density.  Figure  6  and  Table 
1  show  sampling  results  for  two  kinds  of  batteries  having  positive  electrode 
carbon  plates  with  a  thickness  of  0.4  and  0.7  mm.  These  batteries  were 
continuously  and  electrically  discharged  at  an  ambient  temperature  of  20 °C  and 
at  a  current  density  of  10  mA/cm^. 


3.0 


0  CV=3.63  V 


..flectric  discharge 
at  20t.  lOmA/^ 


Voltage 
(V)  2.0 


End  voltage 


2.0  V 


I _ _ _ I _ 1 _ I _ t_ 

0  30  60  90  120  ISO 

Electric  discharge  duration  (  min  ) 


Figure  6.  Electric  Discharge  Characteristics  of 
Table  2  shows  the  electrode  dimension  and  area  of  the  two  batteries. 
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Table  1.  Electric  Discharge  Characteristics  of  Cylindrical  Lithium-Thionyl 
Chloride  Battery 


Thickness  of  positive  electrode  carbon  plate  (mm) 

0.4 

0.7 

Electric  discharge  current  value  (A) 

4.0 

3.5 

Electric  discharge  duration  (hours  and 

minutes) 

2  h  34  m 

2  h  20  m 

Average  voltage  (V) 

3.065 

2.915 

Capacity 

Ah 

10.27 

8.17 

Wh 

31.5 

23.8 

Energy  density 

Wh/kg 

312 

231 

Wh/L 

631 

484 

Table  2.  Dimension  and  Area  of  Electrode 


Thickness  of  positive  electrode 
carbon  plate 

Dimension  and  area  of  electrode 

0 . 4  mm 

500  mm  X  40  mm  =  200  cm^ 

0.7  mm 

430  mm  X  40  mm  =  172  cm^ 

2.2.2  Electric  Discharge  From  Disk-Type  Batteries 

Two  kinds  of  batteries  were  test  manufactured  taking  into  consideration 
resistance  to  vibration  and  impact.  They  were  equipped  with  disk-like 
electrodes.  One  (A  battery)  had  an  effective  electrode  area  of  215  cm^,  while 
the  other  (B  battery)  had  29  cm^.  In  addition,  two  kinds  of  batteries  were 
made  to  prevent  the  generation  of  voltage  delay  phenomena.  The  Li  pole  of  one 
was  coated  with  PVC,  while  that  of  the  other  was  coated  with  CA. 

(1)  Load  current  characteristics 

Figure  7  and  Table  3  show  the  results  of  setting  the  ambient  temperature  of 
the  A  battery  at  20°C  and  continuously  and  electrically  discharging  this 
battery  at  load  currents  of  0.2,  0.5,  1,  2,  and  3  A.  In  the  same  manner, 
Figure  8  and  Table  4  show  the  results  of  continuously  and  electrically 
discharging  the  B  batteries  at  load  currents  of  75,  300,  and  600  mA. 


Ambient 


Figure  7.  Load  Current  Characteristics  of  A  Battery 


Table  3.  Electric  Discharge  Characteristics  of  A  Battery 


Electric  discharge 
current  value  (A) 

Electric  discharge 
duration  (h,  min) 

Average 

voltage 

(V) 

Capacity 

Ah  Wh 

0.2 

26  h  44  min 

3.325 

5.35 

17.79 

0.5 

9  h  16  min 

3.114 

4.63 

14.42 

1 

3  h  35  min 

2.972 

3.58 

10.64 

2 

1  h  23  min 

2.843 

2.77 

7.88 

3 

40  min 

2.575 

2.00 

5.15 

Electric  discharge  (  Ah  ) 


Figure  8.  Load  Current  Characteristics  of  B  Battery 
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Table  4.  Electric  Discharge  Characteristics  of  B  Battery 


Li 

electrode 

treating 

agent 

Electric 
discharge 
current 
value  (mA) 

Electric 
discharge 
duration 
(h,  min) 

Average 

voltage 

(V) 

Capacity 

Ah 

Wh 

CA 

75 

24  h  27  min 

3.250 

1.83 

5.95 

300 

3  h  35  min 

3.178 

1.075 

3.42 

600 

1  h  17  min 

3.112 

0.77 

2.40 

PVC 

75 

18  h  45  min 

3.207 

1.406 

4.51 

300 

2  h  26  min 

3.222 

0.73 

2.35 

600 

48  min 

3.077 

0.48 

1.48 

(2)  Temperature  Characteristics 

Figures  9  and  10  and  Table  5  show  the  results  of  setting  the  ambient 
temperature  of  the  A  battery  at  — 40°C,  — 20°C,  0“C,  20 “C,  40 “C,  and  60 °C  and 
continuously  and  electrically  discharging  these  batteries  at  constant  currents 
of  0.5  and  2  A.  In  the  same  manner,  Figure  11  and  Table  6  show  the  results  of 
testing  the  B  batteries. 


Electric  discharge  duration  (h) 


Figure  9 .  Temperature  Characteris¬ 
tics  of  A  Battery  (Electric 
discharge  at  constant 
current  of  0.5  A) 


I — i — ■...  „i  .  I I  _ >  ■  » 

0  30  60  90  120 

Electric  discharge  duration  (min) 

Figure  10.  Temperature  Characteris¬ 
tics  of  A  Battery  (Electric 
discharge  at  constant 
current  of  2  A) 
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Table  5.  Temperature  Characteristics  of  A  Battery 


Battery  environ¬ 
mental  tempera¬ 
ture  (”0) 

60 

40 

20 

0 

-20 

-40 

Elec¬ 

tric 

dis¬ 

charge 

at 

0.5  A 

Electric 
discharge 
duration 
(h,  min) 

11  h, 

15  min 

10  h, 

12  min 

9  h, 

16  min 

8  h, 

15  min 

3  h, 

25  min 

Average 

voltage 

(V) 

3.379 

3.348 

3.114 

3.007 

2.766 

2.598 

Ca- 

pac 

ity 

Ah 

5.63 

5.10 

4.63 

4.13 

2.57 

1.71 

Wh 

19.0 

17.1 

14.4 

12.4 

7.11 

4.44 

Elec¬ 

tric 

dis¬ 

charge 

at 

2  A 

Elect 
disci 
durat 
(h,  r 

;ric 

large 

:ion 

lin) 

1  h, 

57  min 

1  h. 

50  min 

1  hr, 

23  min 

56  min 

41  min 

13  min 

Average 

voltage 

(V) 

2.937 

2.926 

2.843 

2.620 

2.359 

2.098 

wm 

Ah 

3.90 

3.67 

2.77 

1.87 

1.37 

0.43 

Wh 

11.5 

10.7 

7.88 

4.90 

3.23 

0.90 

Figure  11.  Temperature  Characteristics  of  B  Battery 

(Electric  discharge  at  constant  current  of 
300  mA) 
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Table  6.  Temperature  Characteristics  of  B  Battery  (Electric  discharge 
at  300  mA) 


Battery  environmental 
temperature  (°C) 

60 

40 

-20 

Li 

electrode 

Ca  coating 

Electric  discharge 
duration  (h,  min) 

5  h,  4  min 

3  h,  35  min 

1  h,  37  min 

Average  voltage 
(V) 

3.328 

3.178 

2.921 

Capac 

ity 

Ah 

1.52 

1.07 

0.48 

Wh 

5.06 

3.42 

1.42 

Li 

electrode 

PVC 

coating 

Electric  discharge 
duration  (h,  min) 

4  h,  24  min 

2  h,  26  min 

56  min 

Average  voltage 
(V) 

3.300 

3.222 

2.925 

Capac 

ity 

Ah 

1.32 

0.73 

0.28 

Wh 

4.36 

2.35 

0.85 

(3)  Storage  Characteristics 

Figures  12  and  13  and  Table  7  show  the  results  of  storing  A  batteries  at  room 
temperature  for  6  months  at  a  temperature  of  45°C  for  3  months,  setting  the 
ambient  temperature  at  20°C  and  continuously  and  electrically  discharging  the 
batteries  at  constant  currents  of  0.5  and  2  A. 


Figure  12.  Storage  Characteristics 
of  A  Battery  (Electric 
discharge  at  constant 
current  of  0.5  A) 


0  30  60  90 

Electric  discharge  duration  (min) 

Figure  13.  Storage  Characteristics 
of  A  Battery  (Electric 
discharge  at  constant 
current  of  2  A) 
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Table  7.  Storage  Characteristics  of  A  Battery 


Storage  condition  environmental 
temperature  20 °C 

Immediately 
after  manu¬ 
facturing 

6  months  at 
room  temp¬ 
eratures 

3  months  at 
45'’C 

Electric  discharge 
duration  (h,  min) 

9  h,  16  min 

9  h,  20  min 

7  h,  47  min 

Average  voltage  (V) 

3.114 

3.123 

3.080 

Capacity 

Ah 

4.63 

4.67 

3.89 

Wh 

14.4 

14.6 

12.0 

Elec¬ 
tric 
dis¬ 
charge 
at  2  A 

Electric  discharge 
duration  (h,  min) 

1  h,  23  min 

1  h,  24  min 

1  h,  10  min 

Average  voltage  (V) 

2.843 

2.834 

2.743 

Capacity 

Ah 

2.77 

2.80 

2.33 

Wh 

7.88 

7.94 

6.39 

(4)  Vibration  and  Impact  Characteristics 

Table  9  shows  the  states  (appearance,  shape,  open-circuit  voltage,  etc.)  of 
the  A  battery  obtained  by  impressing  this  battery  under  the  vibration  and 
impact  conditions  shown  in  Table  8.  After  impression,  the  battery  was 
electrically  discharging  at  a  constant  current  of  2  A.  Figure  14  shows  the 
results  of  this  discharge  test. 


Table  8.1  Conditions  of  Vibration 


Division 

I 

II 

Frequency 

10-55  Hz 

55-500  Hz 

Acceleration 

— 

10  G 

Total  amplitude 

2.03  mm 

— 

Sweep  rate 

1  Hz/min 

15  min/cycle 

Direction 

X,Y,Z  (three  axes) 

X,Y,Z  (three  axes) 

Time 

90  rain/axis 

30  min/axis 

Table  8.2  Conditions  of  Impact 


Acceleration 

Time 

Waveform 

Direction 

Number 
of  times 

100  G 

Sine  half  wave 

X,Y,Z 

(three  axes) 

3  times/axis 
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Table  9.  States  Before  and  After  Vibration  and  Impact  Tests 


Division 

Appearance 
and  shape 

Weight 

Open  circuit 
voltage 

Internal 

resistance 

Before 

impression 

No 

abnormality 

108.8  g 

3.65  V 

0.32  n 

After 

impression 

No 

abnormality 

108.8  g 

3.65  V 

0.59  n 

Electric  discharge  duration  (min) 


Figure  lA.  Electric  Discharge  Characteristics  After 
Vibration  and  Impact  Tests 

2.2.3  Safety  of  Aggregated  Batteries 

Two  kinds  (A  and  B)  of  four  disk-type  batteries  were  tentatively  connected  in 
series  to  form  a  group  of  batteries.  The  safety  of  the  batteries  in  this 
group  was  studied  by  subjecting  them  to  the  usual  discharge  tests  together 
with  an  over-discharge  test  and  a  short-circuit  test. 

(1)  Discharge  Characteristics 

Figure  15  and  Table  10  show  the  constant  current  discharge  characteristics  of 
the  group  of  batteries  at  an  ambient  temperature  of  20°C  and  at  a  current 
density  of  about  10  mA/cm^. 


fl  battery 


L . . . . . . . r 

0  12  3  4 

Electric  discharge  duration  (h) 


Figure  15.  Electric  Discharge  Characteristics 
of  Aggregated  Batteries 
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Table  10.  Electric  Discharge  Characteristics  of  Aggregated  Batteries 


A  battery: 
Electric 
discharge  at 

2  A 

Electric  discharge  duration  (h,  min) 

1  h,  55  min 

Average  voltage  (V) 

12.50 

Capacity 

Ah 

3.83 

Wh 

47.9 

B  battery: 
Electric 
discharge  at 
0.5  A 

Electric  discharge  duration  (h,  min) 

3  h,  54  min 

Average  voltage  (V) 

11.77 

Capacity 

Ah 

1.95 

Wh 

23.0 

(2)  Over-Discharge  Characteristics 

Over-discharge  refers  to  a  state  in  which,  even  if  the  terminal  voltage  is 
less  than  the  specified  end  voltage,  the  discharge  will  continue.  A  pole 
change— preventing  diode  was  connected  in  parallel  to  each  unit  battery  to 
prevent  any  abnormalities,  even  if  this  state  continues.  In  addition,  if  poly¬ 
switches  are  connected  in  series  to  a  large-capacity  A  battery,  the  resistance 
of  these  polyswitches  will  increase  in  accordance  with  the  sharp  temperature 
rise  at  the  end  period  of  discharge,  and  the  current  will  be  cut  off.  By  this 
means,  it  was  possible  to  confirm  the  safety  of  the  group  of  batteries.  Also, 
the  polyswitch  has  a  conductive  plastic  element  with  a  temperature  coefficient 
for  which  resistance  is  positive.  No  abnormalities  in  the  appearance  or  shape 
of  this  battery  were  detected. 

Figure  16  shows  the  over-discharge  characteristics  of  an  aggregate  group  of  B 
batteries . 


Electric  discharge  duration  (h) 

Figure  16.  Electric  Over-Discharge  Characteristics 
of  Aggregated  B  Batteries 
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(3)  Short-Circuit  Test  Characteristics 

When  an  A  battery  with  a  large  capacity  shorts  an  external  circuit,  over¬ 
current  will  flow  sharply  and  caloric  power  will  increase.  If  this  state 
continues,  the  safety  valve  of  this  battery  will  activate  itself,  but  in  the 
worst  case,  the  battery  may  explode.  Polyswitches  were  connected  in  series 
with  the  circuit  to  prevent  over-current  from  being  generated  n  this  circuit. 
The  polyswitch  is  a  self-return  device  that  protects  the  element  used  to 
prevent  over-current  and  excessive  temperature  rise  from  being  generated  in 
the  surface  of  each  unit  battery.  Rapid  blowout  fuses  were  connected  in 
series  to  the  B  batteries  and  the  safe  operation  of  these  batteries  was 
checked. 

Figure  17  shows  a  short-circuit  test  circuit  for  the  A  battery.  The  knife 
switch  shown  in  this  figure  was  closed  and  the  battery  voltage,  current,  and 
surface  temperature,  together  with  the  surface  temperature  of  the  polyswitch 
were  measured  with  a  recorder. 
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Figure  17.  Short-Circuit  Test  Circuit 


Figure  18.  Short-Circuit  Test  of  Aggregated  A  Batteries 

Figure  18  shows  the  results  of  a  short-circuit  test  of  the  A  battery.  The 
open-circuit  voltage  and  peak  current  were  3,65  V  and  43  A,  respectively,  but 
immediately  the  voltage  and  current  dropped  to  8  mV  and  about  1.2  A, 
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respectively.  It  was  confirmed  that  the  safety  device  worked  satisfactorily, 
because  the  polyswitch  surface  temperature  and  the  battery  surface  temperature 
were  about  85°C  and  50°C,  respectively. 

The  safety  of  the  B  batteries  was  confirmed  by  the  fact  that  the  fuses  blew 
almost  instantly  during  short-circuit  tests  of  these  batteries. 

3.  Discussion 

3.1  Cylindrical  Battery 

A  feature  of  cylindrical-type  batteries  with  a  spiral  bobbin  electrode  struc¬ 
ture  is  that  it  is  possible  to  widen  the  electrode  area.  Therefore,  these 
batteries  are  most  advantageous  in  terms  of  performance  than  any  other 
batteries.  However,  the  surface  temperature  of  some  batteries  was  60-70°C  at 
the  end  period  of  discharge.  This  is  because  the  radiating  surface  is  smaller 
than  the  discharge  current  value.  Accordingly,  in  order  to  put  such  cylindri¬ 
cal  batteries  to  practical  use  as  high  efficient  discharge  batteries,  their 
safety  must  be  ensured,  given  the  dispersion  structure  of  heat  generated  by  an 
aggregation  of  batteries  and  because  there  is  a  possibility  of  thermorunaway 
being  caused  by  the  thionyl  chloride,  which  has  a  boiling  point  of  77 °C. 

Table  11  shows  a  comparison  of  the  output  densities  for  batteries  with  various 
structures . 


Table  11. 


Comparison  Among  Various  Battery  Structures  in  Respect  of  Output 
Density 

(Current  density,  10  mA/cm^) 


Battery  structure 

Electric 

discharge 

current 

value 

Average 

voltage 

Average 

output 

Output 

density 

W/kg 

W/L 

Cylin¬ 

drical 

battery 

Thickness  of 
carbon  plate 

0.4mm 

4  A 

3.065  V 

12.3  W 

122 

246 

0 . 7mm 

3.5  A 

2.915  V 

10.2  W 

99 

207 

1  Coin  type  battery 

60  mA 

3.133  V 

188  mW 

14 

31 

Disk 

type 

battery 

A  battery 

2.843  V 

5.686  W 

54 

123 

B 

bat¬ 

tery 

Coat¬ 

ing 

CA 

300  mA 

3.178  V 

953. 4mW 

35 

78 

PVC 

300  mA 

3.222  V 

966 . 6mW 

36 

81 
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3 . 2  Disk  Type  Battery 
3.2.1  A  Battery 

(1)  Load  Current  Characteristics 

The  results  depicted  in  Figure  7  and  Figure  3,  show  that  the  flatness  of  the 
discharge  voltage  deteriorated  at  a  current  of  3  A  (current  density:  15  mA/cm^) 
as  compared  with  that  at  ZA  discharge.  The  Wh  capacity  was  65  percent. 
Figure  19  also  shows  the  relationship  between  discharge  current  and  battery 
capacity. 


Figure  19 .  Relationship  Between 
Capacity  and  Electric 
Discharge  Current 

(2)  Temperature  Characteristics 


Figure  20.  Relationship  Between 
Capacity  and  Ambient 
Temperature 


The  results  depicted  in  Figures  9  and  10  and  Table  5  indicate  that  electric 
discharge  at  a  current  of  0.5  A  clearly  affected  the  electric  discharge 
duration  at  temperatures  of  -20°C  or  less,  but  electric  discharge  was  possible 
at  temperatures  as  low  as  — 40°C. 


The  voltage  delay  time  was  more  than  3  minutes  at  an  electric  discharge  of 
2  A  and  at  a  temperature  of  -40°C.  It  was  possible  to  obtain  electric 
discharge  under  these  conditions,  but  in  order  to  achieve  a  highly  efficient 
electric  discharge  at  cryogenic  temperatures,  it  will  be  necessary  to  explore 
ways  to  improve  electrolytes,  coating  agents  for  lithium  pole  surfaces,  etc. 

Figure  20  also  shows  the  relationship  between  ambient  temperature  and  battery 
capacity  (Ah) . 

(3)  Storage  Characteristics 


The  results  depicted  in  Figures  12  and  13  and  Table  7  indicate  that  battery 
capacity  was  not  affected  by  having  been  stored  for  6  months  at  room 
temperature.  But  battery  capacity  after  storage  for  3  months  at  a  temperature 
of  45 °C  was  reduced  by  about  20  percent. 
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(4)  Vibration  and  Impact  Characteristics 

No  abnormalities  were  detected  in  the  results  obtained  from  tests  conducted 
under  the  vibration  and  impact  conditions  shown  in  Tables  8 . 1  and  8.2. 
Further,  as  can  be  seen  from  Table  9,  the  results  obtained  from  subsequent 
electric  discharge  at  2  A  did  not  show  any  effects  from  vibration  and  impact 
within  the  normal  dispersion  of  battery  performance. 

3.2.2  B  Batteries 

Figures  8  and  11  and  Tables  4  and  6  show  the  load  current  characteristics  and 
temperature  characteristics  of  two  kinds  of  batteries  whose  negative  electrode 
Li  pole  surfaces  were  coated  with  PVC  and  CA.  The  results  indicate  that  the 
battery  whose  negative  electrode  Li  pole  surface  was  coated  with  CA  generally 
proved  superior  to  that  coated  with  PVC  in  terms  of  performance. 

3.3  Aggregate  Group  of  Batteries 

Safety  was  checked  by  conducting  over-current  and  short-circuit  tests.  As  a 
result,  it  was  confirmed  that  there  was  neither  any  danger  of  explosion  nor 
any  abnormalities  in  the  batteries,  thanks  to  the  protection  elements 
incorporated  in  an  aggregate  group  of  batteries . 

4.  Conclusion 

Research  on  lithium-thionyl  chloride  batteries  as  primary  batteries  for  highly 
efficient  electric  discharge  has  been  conducted  since  FY  1984.  During  this 
period,  these  kinds  of  batteries  for  electric  discharge  with  a  low  degree  of 
efficiency  have  been  on  the  market  as  power  sources  for  the  backup  of 
integrated  circuit  memories,  etc.,  and  technologies  for  the  safety  of  these 
batteries  have  been  established.  In  the  future,  however,  such  batteries  must 
be  further  studied  and  improved  to  overcome  the  voltage  delay  phenomena 
generated  during  the  initial  stage  of  electric  discharge  at  cryogenic 
temperatures  and  to  solve  the  problem  of  how  to  preserve  highly  efficient  at 
high  temperatures  after  the  batteries  have  been  stored  for  extended  periods. 
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[ Text ]  Summary 

We  manufactured  a  spherical  lens  antenna  using  a  circular  planar  spiral 
antenna  covering  8~18  GHz  range  and  measured  its  characteristics. 

As  the  result  of  our  measurement,  the  antenna  gain  and  beam  width  were 
determined  to  be  about  14-24  dBi  and  5-12  degrees  in  a  range  from  8-18  GHz. 
The  angular  error  of  arrival  was  0 . 3-1 . 2'’rms ,  based  on  model  analysis  using 
the  measured  antenna  pattern. 

The  results  of  our  measurement  and  model  analysis  are  expected  to  be  applied 
to  the  design  of  a  lens  antenna  using  electronic  support  measures. 

1.  Introduction 

This  report  describes  the  contents  of  what  was  primarily  an  experimental  study 
of  a  lens  antenna  equipped  with  a  receiving  dielectric  sphere  (hereafter 
called  the  "lens").  This  antenna  has  been  used  in  two-dimensional  angle 
measuring  work  as  part  of  the  "Research  on  ESM"  program  conducted  by  TRDI's 
Second  Research  Center  in  FY  1988. 

There  is  an  extensive  body  of  research  on  using  lens  antennas  as  optical 
scanning  antennas.  A  lens  antenna  is  designed  so  that  a  focal  point  is 
established  on  the  sphere  or  its  external  portion  by  using  a  spherical 
refracting  substance.  Therefore,  wave  fronts  incident  from  arbitrary 
directions  converge  at  this  focal  point. ^ 
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Such  antennas  possess  the  following  features; 

(1)  A  niomber  of  equivalent  phases  can  be  formed  in  a  direction 
symmetrical  to  the  lens  sphere  by  arranging  a  number  of  element  antennas 
on  or  along  the  lens  sphere ;  these  can  be  used  as  a  number  of  beam 
antennas . 

(2)  When  a  number  of  element  antennas  are  arranged  on  the  surface  of 
the  sphere  or  along  its  external  portion,  the  converging  area  will 
virtually  become  two-dimensional  hemisphere  as  long  as  incoming  waves 
are  not  blocked  by  element  antennas  arranged  on  the  opposite  side  of  the 
sphere . 

(3)  The  frequency  band  is  wide;  this  is  because,  in  general,  it  is 
possible  to  use  a  dielectric  sphere  where  the  change  in  the  dielectric 
constant  is  less  than  that  of  the  frequency. 

When  these  features  are  incorporated  into  this  antenna,  which  has  a  simple 
structure,  the  antenna  appears  to  be  promising  as  a  wide  band  receiving 
antenna  for  two-dimensional  angle  measuring  work.^ 

Previously,  the  author  has  theoretically  studied  the  wide  band  properties  and 
angle  measuring  properties  of  the  antenna.^  In  this  research  project,  a  lens 
antenna  was  made  by  installing  a  spiral  antenna  (8-18  GHz)  as  a  wide  band 
element  antenna  on  a  lens  with  a  diameter  of  177  mm.  The  wide  band  properties 
and  angle  measuring  properties  of  this  lens  antenna  were  studied 
experimentally,  and  the  angle  measuring  accuracy  was  simulated  by  using  a 
result  obtained  from  this  study. 

As  a  result,  it  was  confirmed  that  this  lens  antenna  possessed  a  wide  band 
property  with  a  gain  of  14—24  dBi  in  the  X  band  and  Ku  band  (8—18  GHz)  .  It  was 
also  confirmed  experimentally  that  the  half  beam  width  of  the  antenna  pattern 
was  as  sharp  as  5  to  about  12  degrees  and  the  first  side  lobe  was  about  17  dB 
as  compared  with  the  maximum  gain  of  the  main  lobe.  In  addition,  it  was 
determined  that  it  is  possible  to  carry  out  angle  measuring  work  with  a  high 
degree  of  accuracy,  because  when  the  angle  measuring  accuracy  was  simulated  by 
using  the  above  result,  the  angle  measuring  error  was  0.3-1.2°rms.  We  will 
hereunder  describe  the  details  of  our  research  project. 

2 .  Antenna 

Photographs  1  and  2  [not  reproduced]  show  the  appearance  of  the  lens  antenna 
made  for  this  research  project  and  that  of  the  spiral  antenna  used  therein. 
As  shown  in  these  photographs,  the  lens  used  in  our  research  project  had  a 
diameter  of  177  mm.  Its  refractive  index  is  distributed  almost  in  accordance 
with  Luneberg's  solution,  shown  as  a  broken  line  in  Figure  1,  by  setting  a 
radius  standardized  from  the  center  of  the  sphere  on  the  transverse  axis  shown 
in  the  figure  and  setting  the  refractive  index  on  the  ordinate  shown  in  the 
figure,  and  by  approximately  constituting  the  lens  with  seven  layers  as  an 
indication  of  the  distribution  of  the  refractive  index. 4,  Also,  as  can  be  seen 
from  this  distribution  of  the  refractive  index,  the  focal  point  is  constructed 
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Figure  1.  Distribution  of  Refractive  Index  of  Lens 

so  that  it  is  positioned  almost  on  the  surface  of  the  sphere.  The  materials 
were  made  by  mixing  Styrofoam  with  very  thin  aluminum  wires  to  control  the 
refractive  index  and  to  minimize  loss . 

The  antenna  was  constituted  by  installing  a  wide  band  spiral  antenna  with  a 
diameter  of  20  mm  as  an  element  antenna  (shown  in  Photograph  2  [not 
reproduced] )  on  this  lens  so  that  it  is  possible  continuously  to  move  this 
spiral  antenna  from  the  surface  of  the  sphere  in  the  direction  of  the  radius. 
Appendix  A-1  shows  the  results  of  measuring  the  gain  of  the  element  antenna 
used  in  the  research  project.  These  results  indicate  that  the  band  of  the 
element  antenna  almost  covers  a  range  of  8-18  GHz. 

Also,  a  lens  installing  table  and  an  element  antenna  installing  table  were 
made  by  using  Styrofoam,  and  thus  the  influence  of  electromagnetic  fields  on 
these  tables  was  minimized. 

3 .  Measuring  Method 

Figure  2  shows  a  system  for  measuring  gain  and  antenna  patterns.  When  they 
were  actually  measured  using  this  system,  as  shown  in  this  figure,  vertically 
and  horizontally  polarized  waves  were  irradiated  from  a  transmitting  antenna 
in  an  electric  dark  room  and  these  waves  were  received  through  a  receiver  by 
means  of  the  standard  antenna  and  the  main  antenna  installed  on  a  rotation 
table.  The  results  were  indicated  with  a  pattern  recorder  and  the  absolute 
gain  as  found  by  comparing  this  receiving  level  with  the  known  gain  of  the 
standard  horn.  Also,  the  element  antenna  was  used  to  receive  circularly 
polarized  waves.  There  is  a  possibility  that  a  polarization  dependence  could 
exist  in  the  element  antenna  itself.  The  clamp  face  of  this  element  antenna 
was  fixed  at  all  times  and,  as  mentioned  above,  the  gain  for  each  respective 
case  was  measured  and  compared  with  every  other.  This  was  done  by  changing  the 
polarized  waves  of  the  transmitting  antenna  in  order  to  check  this  particular 
characteristic . 
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Figure  2.  Block  Diagram  of  Measuring  System 


4.  Measurement  Results  and  Discussion 

The  results  of  measuring  the  gain  and  antenna  patterns  of  the  antenna  made  for 
this  research  project  are  shown  below.  Also,  the  following  measurement  results 
were  obtained  by  closely  attaching  the  element  antenna  to  the  surface  of  the 
lens . 

4.1  Antenna  Patterns 

Figure  3(a)  to  (c)  is  an  example  of  the  results  of  measuring  the  antenna 
patterns  of  the  horizontally  polarized  waves .  These  patterns  were  measured  by 
setting  the  elevation  angle  to  0®  and  by  changing  the  azimuth  angle.  The 
transverse  and  the  ordinate  denote  the  azimuth  angle  and  the  gain  (dBi) , 
respectively,  in  (a)  to  (c) . 

As  a  result,  at  8-18  GHz,  the  half  beam  width  ranged  from  5-12°.  Sharp  beams 
extending  over  a  wide  band,  as  compared  with  the  maximum  gain  of  the  main  lobe 
were  obtained;  that  of  the  first  side  lobe  was  restricted  to  about  17  dB,  and 
wide  band  properties  were  confirmed. 

As  mentioned  above,  a  sharp  beam  width  was  obtained  at  each  frequency,  but  in 
order  to  obtain  sharper  beam  widths,  it  seems  to  be  necessary  to  theoretically 
determine  the  distribution  of  phases  and  the  amplitude  of  antenna  patterns  of 
the  element  antenna  in  reverse  form  the  optimum  distribution  of  electro¬ 
magnetic  fields  of  the  aperture  of  the  lens  antenna  and  to  use  element 
antennas  possessing  such  characteristics. 

4.2  Gain 

Figure  4  shows  the  results  of  measuring  the  gains  of  horizontally  and 
vertically  polarized  waves .  The  transverse  and  the  ordinate  denote  the 
frequency  (GHz)  and  the  gain  (dBi),  respectively,  in  this  figure. 

As  a  result,  it  was  ascertained  that  a  gain  of  14-24  dBi  could  be  obtained  at 
frequencies  of  8-18  GHz.  Also,  as  previously  mentioned,  it  can  be  expected 
that  if  the  gain  of  the  element  antenna  and  the  distribution  phases  are 
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figure,  a  comparison  of  the  results  obtained  from  this  measurement  indicates 
that  the  difference  is  about  0.5  dB  at  most  and  that  there  is  little 
dependence  of  the  waves  on  the  gain. 

5.  Simulation  of  Angle  Measuring  Accuracy 

The  angle  measuring  accuracy  of  this  lens  antenna  was  studied  as  an  example  of 
an  angle  measuring  method  based  on  an  assumption  concerning  an  amplitude 
comparison  system  employing  a  direct  wave  detecting  system  that  will  be 
mentioned  later  on. 

In  general,  when  the  angle  measuring  accuracy  of  a  system  is  evaluated,  it  is 
necessary  to  study  receiving  systems  and  various  methods  of  correcting  angle 
measuring  errors  as  well  as  antenna  systems.  In  this  case,  the  following  can 
be  cited  as  principal  error  factors  that  affect  angle  measuring  accuracy: 

(a)  Change  of  characteristics,  accompanying  temperature  change,  etc., 
of  the  various  components  that  constitute  the  receiver. 

(b)  Fluctuation  of  gain  between  a  pair  of  antennas. 

(c)  Change  of  antenna  patterns  accompanying  frequency  change. 

(d)  Change  of  antenna  patterns  accompanying  change  of  polarized  waves. 

Taking  into  account  these  error  factors,  the  following  hypotheses  were  studied 
as  examples  of  simulations  related  to  angle  measuring  accuracy  when  the  angle 
measurements  are  carried  out  using  this  antenna. 

(i)  It  shall  be  possible  to  memorize  error  correction  factors 
accompanying  items  (b)  and  (c)  above,  which  are  error  factors  of  the 
antenna  system,  for  every  frequency  in  advance  and  to  correct  these 
errors  during  the  measurement  of  angles . 

(ii)  As  previously  mentioned,  the  angle  measuring  system  shall  be  an 
amplitude  comparison  system  employing  a  direct  wave  detecting  receiver 
(a  block  diagram  of  this  system  is  shown  in  Appendix  A-2) . 

(iii)  Element  antennas  shall  be  arranged  on  the  respective  antennas  so 
that  these  antennas  cross  each  other  at  a  power  gain  mesial  point 
decreased  by  three  decibels. 

As  hypothesized  in  item  (i),  when  considering  that  errors  caused  by  the 
antenna  system  are  corrected  during  the  measurement  of  angles,  the  error 
factor  to  be  considered  can  be  limited  only  to  item  (a)  above.  Then,  assuming 
that  the  rms  value  of  each  component  error  of  these  receivers  is  Aa,- ,  rms 
value,  AA  of  the  error  of  all  the  receivers  is  given  in  the  following 
equation: 


AA  = 


(1) 
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Next,  calculating  the  difference  in  gain  between  a  pair  of  antennas,  drawing 
an  S  curve  on  the  basis  of  the  result  of  this  calculation,  and  finding  the 
inclination  of  a  of  this  curve,  the  angle  measuring  error  is  found  by  using 
AA  obtained  from  equation  (1)  and  is  expressed  as  the  following  equation: 

A0  =  AA/a  (2) 

First  the  S  curve  was  determined  by  measuring  the  antenna  patterns  shown  in 
Figure  3  in  view  of  the  above  theories  in  order  to  simulate  the  angle 
measuring  error.  Figure  5  shows  an  example  of  the  S  curve.  Finding  the 
inclination  a  neighboring  0“  from  these  figures  gives  Table  1. 


II 


Figure  5 .  Example  of  S  Curve 


Table  1.  Inclination  of  S  Curve  Determined  by  Measured  Antenna  Patterns 


Frequency  (GHz) 

8 

14 

18 

Inclination  of  S 
curve  (dB/deg) 

1.2 

4.0 

5.8 

Also,  errors  generated  by  the  general  wave  detector,  logarithmic  amplifier, 
etc.,  are  as  shown  in  Table  2  on  the  basis  of  references  5,  6,  etc. 

When  the  angle  measuring  error.  Ad,  is  calculated  by  using  equations  (1)  and 
(2)  taking  into  consideration  the  above  results,  the  results  shown  in  Table  3 
will  be  obtained. 
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Table  2.  Dispersion  of  Characteristics  of  Single  Substance  Receiver  and 
Between  Receivers 


Error 

factor 

Standard  deviation  (dB) 

Wave  detector 

Tracking 

0.7 

Linearity 

0.7 

Logarithmic  amplifier 

Tracking 

0.7 

Linearity 

0.7 

Cables 

Tracking 

0.35 

Total 

1.44 

Table  3.  Result  of  Calculating  Angle  Measurement  Errors 


Frequency  (GHz) 

8 

14 

18 

Angle  measuring 
error  (°  rms) 

1.2 

0.4 

0.3 

As  previously  mentioned,  the  above  hypotheses  were  studied  as  examples  of 
simulations  related  to  angle  measuring  accuracy.  It  is  believed  that  even  if 
one  receiving  system  is  changed  to  another,  an  angle  measuring  accuracy  close 
to  the  values  obtained  from  the  simulation  can  be  obtained. 

6.  Conclusion 

In  this  research  project,  we  created  a  receiving  lens  antenna  that  was  used  to 
measure  two-dimensional  angles.  This  antenna  was  studied  by  conducting 
experiments  and  simulations  to  determine  the  wide  band  properties  and  angle 
measurement  properties  of  this  antenna. 

As  a  result  of  this  study,  it  was  confirmed  that  this  antenna  had  wide  band 
properties  with  a  gain  of  14-24  dBi  in  the  X  and  Ku  bands  (8-18  GHz).  It  was 
further  confirmed  experimentally  that  the  half  beam  width  of  the  antenna  beams 
was  as  sharp  as  5  to  about  12°  and  that  the  maximum  gain  of  the  first  side 
lobe  was  about  17  dB  lower  than  that  of  the  main  lobe .  As  a  result  of  conduct¬ 
ing  experiments  and  simulations  on  angle  measuring  accuracy  on  the  basis  of 
this  result,  it  was  ascertained  that  it  is  possible  to  measure  angles  with  a 
high  degree  of  accuracy  (the  angle  measurement  error  is  0.3-1. 2°  rms) . 

These  results  would  appear  to  be  useful  for  the  manufacture  of  practical 
receiving  lens  antenna  to  be  used  to  measure  two-dimensional  angles. 
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Appendix  A.l  Gain  of  Element  Antenna 


Figure  A.l  shows  the  results  of  measuring  gain  against  the  horizontally 
polarized  waves  of  element  antennas  used  in  this  research  work. 


Figure  A.l  Results  of  Measuring  Gain  of  Element  Antennas 
(Horizontally  polarized  waves) 

Appendix  A. 2  Receiving  System 

Figure  A. 2  shows  a  block  diagram  of  a  direct  wave  detecting  and  receiving 
system  hypothesized  as  a  receiving  system. 


Antenna 


Figure  A. 2  Block  Diagram  of  Direct  Wave  Detection  and 
Receiving  System 
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Motor  (I) — Effectiveness  of  Two-Stage  Thrust  System  and  Characteristics  of 
Porous  Bulkhead  Mechanism" ] 

[ Text ]  Summary 

Combustion  tests  were  conducted  in  order  to  develop  an  optimized  two-stage 
rocket  motor.  The  pulse  motor  consisted  of  two  propellant  grains,  two 
igniters,  and  a  convergent-divergent  nozzle.  Both  chambers  were  separated  with 
a  disk  which  was  partially  opened  when  the  sustainer  chamber  was  operated. 
Test  results  Indicated  that  a  sustainer  pulse  was  generated  at  a  designated 
time  after  a  booster  pulse  operation.  The  two-pulse  technique  developed  in 
this  study  was  found  to  be  very  useful  for  practical  application. 

1.  Introduction 

The  interception  capabilities  of  missiles  can  be  expanded  by  improving  the 
performance  of  their  solid  rocket  motors.  Up  to  now,  efforts  to  enhance  solid 
rocket  motor  performance  have  focused  on  increasing  the  specific  impulse  of 
propellants  and  lightening  motor  cases.  However,  we  have  reached  the  point 
where  it  is  difficult  to  increase  performance  significantly  by  using  only 
these  technologies.  Meanwhile,  lift-fire  technologies  for  antiaircraft 
missiles  have  responded  enthusiastically  to  the  challenges  posed  by  the  rapid 
rise  in  maneuverability,  low  altitude  cruising  capability,  etc.,  of  aircraft. 
In  order  to  meet  such  requirements,  it  is  necessary  to  raise  the  overall 
performance  of  rocket  motors.  Accordingly,  the  thrust-control  rocket  motor, 
which  improves  performance  by  controlling  thrust,  can  be  regarded  as  a  new 
system  to  be  used  in  conjunction  with  other  high-performance  technologies.^ 
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This  system  is  designed  so  that  air  resistance  proportional  to  second  power  of 
velocity  can  be  reduced  by  restraining  the  initial  velocity  of  the  missile; 
lift-fire  and  terminal  velocity  can  be  increased  and  maneuverability  can  be 
ensured’^  by  controlling  the  operation  of  the  rocket  motors.  In  this  research 
work,  a  small  thrust-control  rocket  motor  was  trial  manufactured,  and  a 
destructive  bulkhead  system,  which  is  a  main  technology,  was  established  by 
conducting  combustion  tests  and  analyses. 

2.  Outline  and  Structure  of  Thrust -Control  Rocket  Motor 

2.1  Outline 

Japanese  defensive  missiles  are  subjected  to  atmospheric  drag  at  all  times 
during  their  flight,  because  they  operate  in  the  atmosphere. 

The  rocket  motors  used  in  the  present  generation  of  surface-to-air  missiles 
[SAMs]  and  air-to-air  missiles  [AAMs]  are  designed  to  use  the  residual 
velocity  when  these  missiles  are  used  to  intercept  the  enemy  by  enhancing  the 
initial  velocity.  However,  the  use  of  this  method  creates  severe  air 
resistance  proportional  to  the  second  power  of  velocity,  and  results  in  large 
energy  losses,  because  of  the  relationship  between  air  resistance  (D)  and 
flight  velocity  (V) ,  as  shown  in  the  following  equation: 

D  =  (1) 


Where , 


Cp:  Coefficient  of  drag 
A  :  Project  area  of  front  of  missile 
p  :  Air  density 

By  contrast,  when  the  initial  velocity  is  lowered,  it  becomes  difficult  to 
ensure  that  sufficient  velocity  remains  to  achieve  interception.  It  is 
believed  that  thrust  can  be  adequately  controlled  by  using  a  rocket  in  which 
it  is  possible  to  control  the  thrust  continuously.^  Thanks  to  this  control, 
it  will  become  possible  to  make  missiles  fly  by  restraining  their  initial 
velocities  and,  by  lowering  the  loss  of  energy  during  their  initial  flights, 
to  enable  them  to  fly  at  their  residual  velocities  during  their  intermediate 
flight  and  to  ensure  interception  velocities  by  raising  the  thrust  and 
reaccelerating  them  during  interception.  This  means  that  it  will  become 
possible  to  increase  lift-fire  and  terminal  velocity,  to  ensure^  a  high  degree 
of  maneuverability,  and  to  expand  the  missile's  interception  range. 

The  liquid-fuel  rocket  engine  can  be  cited  as  a  rocket  in  which  it  is  possible 
to  control  thrust  continuously.  This  rocket  employs  a  liquid  fuel  as  a 
propellant  and  it  is  possible  to  control  the  thrust  by  controlling  the  amount 
of  propellant  burned  in  the  rocket.  In  addition,  the  maximum  specific  impulse 
(Isp)  is  410  seconds,^  which  is  higher  than  that  (about  260  seconds)  of  solid 
rocket  motors.  Prompt  firing  properties,  handling  properties,  and  ease  of 
storage  are  indispensable  elements  for  tactical  missiles.  Liquid-fueled 
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rockets  present  problems  in  this  regard,  because  the  rocket  has  a  complex 
structure^  and  it  is  difficult  to  handle  the  liquid  fuel.  By  contrast,  solid- 
fuel  rockets  excel  in  this  respect,  and  are  used  in  most  of  the  missiles 
currently  in  the  inventories  of  the  advanced  Western  countries.  At  present, 
solid— fuel  missiles  are  gradually  replacing  liquid— fueled  ones. 

The  disadvantage  of  the  solid-fuel  missile  is  that  it  is  difficult  to  variably 
control  the  thrust,  because  it  is  impossible  to  shut  off  the  solid  rocket 
motor  during  the  combustion  of  the  solid  propellant  or  to  control  the 
combustion  rate  of  the  solid  propellant.®  However,  use  of  the  following 
systems  may  make  it  possible  to  achieve  thrust-control  rocket  motors  suitable 
for  defensive  tactics,  because  these  systems  control  thrust  by  dividing  the 
propellant . 

When  solid  rocket  motors  are  used  in  missiles,  it  should  be  possible  to 
control  the  thrust  hy  dividing  the  propellant  into  some  nvimber  of  parts, 
controlling  the  respective  ignition  times,  intermittently  controlling  the 
operation  and  optimizing  the  timing  of  the  operation.  The  multistage  (two- 
stage)  rocket  shown  in  Figure  1(a)  and  the  multistage  (two-stage)  thrust 
rocket  shown  in  Figure  1(b)  can  be  regarded  as  examples  of  such  systems. 
Independent  rocket  motors  are  connected  in  series  in  the  multistage  rocket, 
and  the  inside  of  the  rocket  motor  is  divided  into  several  combustion  chambers 
in  the  multistage  thrust  rocket. 


1  Second  stage  S&A  device 

2  Second  stage  igniter 

3  Second  stage  propellant 

4  Second  stage  nozzle 

5  Separator 

6  First  stage  S&A  device 

7  First  stage  igniter 

8  First  stage  propellant 

9  First  stage  nozzle 


1  Second  stage  S&A  device 

2  Second  stage  igniter 

3  Second  stage  propellant 

4  Bulkhead 

5  First  stage  igniter 

6  First  stage  propellant 

7  Nozzle 


Figure  1(a).  Two-Stage  Rocket  Figure  1(b).  Two-Stage  Thrust  Rocket 

The  use  of  the  former  system  will  improve  propulsion  efficiency,  because 
unnecessary  mass  is  jettisoned  by  the  orderly  separation  of  spent  rockets  from 
the  missile.^  Further,  it  is  possible  to  enhance  propulsion  efficiency  by 
optimizing  the  nozzle  shape  corresponding  to  the  pressure  of  each  combustion 
chamber,  because  a  special-purpose  nozzle  is  installed  on  each  rocket  motor. 
However,  this  system  is  inappropriate  for  tactical  missiles,  for  the  following 
two  reasons:  1)  the  possibility  that  the  rockets  separated  from  the  missile 
will  fall  to  the  ground,  and  2)  a  number  of  nozzles  with  large  capacities  are 
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installed  in  the  rocket  motors,  while  the  amount  of  propellant  is  sharply 
limited  because  the  outside  dimensions  of  these  missiles  are  severely 
restricted.  Therefore,  this  system  is  adopted  for  space  rockets  whose  launch 
sites  are  determined  in  advance  and  whose  capacities  can  be  ensured. 

The  latter  system  has  no  separating  section  and  the  structure  of  this  system 
is  simplified.  The  system  is  designed  so  that  a  number  of  combustion  chambers 
use  a  common  nozzle,  and  the  combustion  gas  generated  in  the  combustion 
chambers  at  the  rear  stage  passes  through  the  inside  of  the  spent  combustion 
chambers.  It  is  also  possible  to  prevent  a  decrease  in  the  rate  of  propellant 
charged  for  the  capacity  of  the  rockets.  A  large  amount  of  combustion  gas  is 
generated  at  high  temperatures  in  the  combustion  chamber  operated  previously. 
When  the  intervals  between  the  respective  combustion  chambers  are  considered, 
it  is  necessary  to  isolate  them  from  each  other  so  that  the  heat,  pressure, 
and  flames  caused  by  the  combustion  gas  passing  through  them  do  not  ignite 
unburned  propellant  stored  in  unoperated  combustion  chambers.  In  addition, 
when  the  combustion  chambers  at  the  rear  stage  are  operated,  the  system  must 
be  structured  so  that  the  combustion  chambers  at  the  nozzle  side  will  no 
longer  be  isolated.  Accordingly,  it  is  necessary  to  solve  a  number  of 
technical  problems  related  to  the  isolation  structure  of  the  bulkhead,  etc., 
which  will  satisfy  these  conditions.  However,  when  this  system  is  compared  to 
the  former  system,  the  latter  system  is  more  suitable  as  a  rocket  motor  for 
tactical  missiles  than  the  former. 

For  these  reasons ,  thrust  control  is  considered  in  this  research  work  in  terms 
of  a  system  in  which  the  inside  of  a  solid  rocket  motor  is  divided  into 
several  combustion  chambers. 

2 . 2  Comparison  With  Current  Rocket  Motors 

Figure  2  shows  a  comparison  between  typical  rocket  motors  used  in  operational 
missiles  and  thrust-control  rocket  motors  with  regard  to  structures  and  thrust 
patterns . 

Figure  2(a)  shows  the  structure  and  thrust  curve  of  the  rocket  motor  used  in 
operational  missiles.  The  rocket  motor  consists  of  a  nozzle  and  a  double 
rear-shaped  propellant  formed  by  overlapping  two  kinds  of  propellants .  The 
relation  between  thrust  and  time  is  that  a  large  thrust  is  obtained  from  the 
combustion  of  the  booster  propellant  on  the  inside  while  a  lesser  thrust  is 
generated  from  the  combustion  of  sustainer  propellant  on  the  outside. 

Figure  2(b)  shows  the  case  of  a  thrust-control  rocket  motor.  This  rocket 
motor,  which  has  a  two-stage  thrust  structure,  consists  of  a  nozzle  and  two 
combustion  chambers.  Different  propellants  can  be  loaded  in  the  respective 
combustion  chambers.  The  two  combustion  chambers  and  the  propellant  are 
separated  from  each  other  with  a  bulkhead  that  shuts  out  the  heat,  flames,  and 
pressure  of  the  combustion  gas.  The  rocket  motor  operating  system  is  designed 
so  that  after  the  first  stage  propellant  is  burned  and  the  arbitrarily  set 
time  passes,  the  second  stage  propellant  will  ignite  and,  simultaneously, 
second  stage  thrust  will  be  generated  by  breaking  the  pressure  bulkhead. 
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Figure  2(a).  Thrust  Curve  and  Figure  2(b).  Thrust  Curve  and 

Structure  of  Operational  Structure  of  Thrust 

Rocket  Motor  Control  Rocket  Motor 


Therefore,  the  thrust  generated  from  this  rocket  motor  is  characterized  by  the 
fact  that  this  thrust  has  two  kinds  of  thrust  depending  on  time,  and  the  time 
interval  between  the  first  and  second  stages  can  be  determined  arbitrarily. 


Figure  3.  Flight  Performance  Using  Thrust-Control 
Rocket  Motor 
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Figure  3  shows  the  relationship  between  flight  velocity  and  distance  for  SAMs 
when  operational  rocket  motors  and  thrust-control  rocket  motors  are  used  in 
SAMs,  respectively.  When  the  operational  rocket  motor  is  used  in  a  SAM,  the 
flight  velocity  will  reach  its  peak  at  the  time  when  the  propellants  of  both 
the  booster  and  sustainer  are  burned  completely,  and,  subsequently,  the  SAM 
will  glide  freely  until  it  hits  the  target.  The  high  maximum  velocity  means 
that,  as  shown  in  equation  (1),  a  large  amount  of  kinetic  energy  is  lost 
because  the  air  resistance  imposed  on  the  rocket  motor  operating  system  is 
proportional  to  the  square  of  the  velocity.  By  contrast,  when  a  SAM  is 
launched,  a  thrust-control  rocket  motor  will  burn  only  the  first  stage 
propellant,  and  thereby  restrain  its  initial  maximum  velocity.  Subsequently, 
this  rocket  motor  will  interrupt  its  operation  by  itself  and  the  SAM  will 
glide  freely.  In  addition,  when  the  velocity  of  the  SAM  reaches  minimum 
flight  velocity,  the  rocket  motor  will  ignite  the  second  stage  propellant  and 
will  again  enhance  flight  velocity.  After  combustion  of  the  second  stage 
propellant  is  completed,  the  SAM  will  glide  freely  until  it  hits  the  target. 
As  mentioned  above,  air  resistance  can  be  reduced  by  restraining  the  maximum 
velocity,  and  lift  fire,  as  shown  in  Figure  3,  can  be  effected  qualitatively 
by  minimizing  energy  loss.  Also,  when  a  SAM  approaches  its  target,  and  when 
its  residual  velocity  is  comparatively  high,  the  ignition  of  the  second  stage 
rocket  motor  can  be  initiated  arbitrarily,  and  it  will  thus  be  possible  to 
increase  the  intercept  velocity  of  the  SAM  by  reaccelerating  it  at  an  early 
period  and  to  ensure  its  maneuverability. 
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Figure  4(a).  Flight  Curve  of  Double 
Rare  Propellant  Rocket 
Motor 


Figure  4(b).  Flight  Curve  of  Two- 
Stage  Thrust-Control 
Rocket  Motor 


Figure  4  shows  the  results  obtained  from  theoretical  calculations  of  the 
relationship  between  the  ballistic  flying  range  and  altitude  of  a  SAM. 
(Figure  4(c)  schematically  depicts  the  results  obtained  from  simulation 
calculations.)  When  a  rocket  motor  (conventional  dual  thrust  motor)  employing 
a  double  rare  propellant  is  used  in  a  missile  on  the  assumption  that  the 
interception  parameters  are  at  an  altitude  of  4,000  or  more  and  the  velocity 
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•  _  ,  T 

U  —  — qw  +  - Cx  —  g-sin0  +  — 

■'  II 

w  —  qu  +  — — C2  4-  gcos0 


X  =  U  COS0  +  W  sin0 


Z  —  u  sin0  —  w  COS0 

Where , 

u,  w:  Velocity  (body  axis) 

q:  Pitch  rate  (horizontal  flying  transfer  90‘’/s) 

Cx.Cj:  Aerodynamic  coefficient  (application  of  aerodynamic 
characteristics  of  usual  SAM,  for  example,  Cxo  =0.6 
at  M  =  2) 

Q:  Dynamic  pressure 
S:  Standard  area  (0.071  m^) 

m:  Mass  (initial  period,  525  kg/g;  propellant,  225  kg/g) 
T:  Thrust  (5,000  kgf  x  6  s-interruption  15  s~5,000  kgf 
X  5  s) 

9 :  Airframe  attitude  angle 
X,Z:  Inertial  coordinate  system 
g:  Gravitational  acceleration 

Figure  4(c) .  Basic  Equation  of  Flight  Simulation 
and  Example  of  Input  Data 

is  400  or  more  per  second,  as  shown  in  Figure  4(a) ,  the  maximum  range  will  be 
about  35  km.  When  a  thrust-control  rocket  motor  is  used  instead  of  a  dual 
thrust  motor  in  this  missile,  and  when  the  gliding  flight  of  the  missile  is 
calculated,  as  shown  in  Figure  4(b),  the  maximum  range  will  be  about  43  km. 
In  this  case,  the  dimensions  and  amount  of  propellant  for  this  thrust-control 
rocket  motor  must  be  equivalent  to  those  of  the  former  rocket  motor.  This 
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fact  indicates  that  it  is  possible  to  increase  the  lift  fire  by  about  25 
percent. 

As  mentioned  above,  it  is  possible  to  raise  the  overall  performance,  including 
lift  fire  and  intercept  capability,  of  missiles  by  using  thrust-control  rocket 
motors  in  these  missiles. 

2.3  Structure  and  Operating  System 

Figure  5  shows  the  structure  of  a  two-stage  thrust-control  rocket  motor  that 
consists  primarily  of  two  combustion  chambers  for  thrust  control,  a  bulkhead 
to  serve  as  a  boundary  between  the  two  combustion  chambers,  a  nozzle,  and  a 
rocket  operating  igniter.  It  is  possible  to  load  different  propellants  in  the 
respective  combustion  chambers. 


1  Second  stage  S&A  device 

2  Second  stage  igniter 

3  Second  stage  propellant 

4  Bulkhead 

5  First  stage  S&A  device 

6  First  stage  igniter 

7  First  stage  propellant 

8  Nozzle 

Figure  5.  Structural  Drawing  of  Two-Stage 
Thrust-Control  Rocket  Motor 

The  bulkhead  must  be  structured  so  that  heat,  pressure,  and  flames  are  iso¬ 
lated  from  each  other  and  so  that  the  second  stage  propellant  is  not  ignored 
by  the  high-temperature  gas  generated  in  the  first  stage  combustion  chambers 
during  the  operation  of  this  combustion  chamber.  It  must  also  be  structured 
so  that  it  can  be  broken  readily  by  second  stage  combustion  and  so  that  the 
nozzle  will  not  become  clogged  with  the  fragments  created  when  it  breaks. 

It  is  necessary  to  consider  lightening  the  safety  and  arming  [S&A]  mechanism, 
which  consists  of  an  igniter  and  a  safety  device,  and  the  position  for 
installing  this  mechanism,  because  the  igniter  and  the  safety  device  have 
independent  combustion  chambers  join  the  S6tA  mechanism.  In  particular,  the 
igniter  for  the  first  stage  combustion  chamber  may  be  regarded  as  a  burning- 
type  igniter  after  the  work  is  finished  so  that  it  will  not  cause  any  serious 
trouble,  such  as  clogging  of  the  nozzle,  etc.,  during  the  operation  of  the 
second  stage  combustion  chamber. 
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Ignition  of  first 
1  stage  propellant 


3  Interruption  of  combustion 


4  ‘i 


Ignition^of  second 


5  Breakage  of  bulkhead 


Figure  6.  Operating  System  of  Thrust-Control  Rocket  Motor 

As  shown  in  Figure  6,  the  second  stage  operating  system  of  the  thrust-control 
rocket  motor  conforms  to  the  following  procedure: 

(1)  The  first  stage  propellant  igniter  fires. 

(2)  Initial  thrust  is  generated  by  burning  the  first  stage  propellant. 
At  this  time ,  the  bulkhead  shuts  out  combustion  gas  and  heat  from  the 
other  chamber. 

(3)  Next,  operation  of  the  rocket  motor  is  stopped  for  an  arbitrary 
time  lasting  from  a  few  seconds  to  10  seconds  or  more. 

(4)  The  second  stage  propellant  igniter  is  operated  in  accordance  with 
instructions  given  from  the  ground  or  from  the  missile  itself. 

(5)  The  second  stage  propellant  is  burned.  The  bulkhead  is  broken  by 
the  combustion  chamber  pressure  generated  at  this  time. 

(6)  Secondary  thrust  is  generated. 
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3.  Method  of  Experiment  and  Device 

The  thrust-control  rocket  motor  involves  three  technical  problems:  a  bulkhead 
breaking  system,  pressure  resistance,  and  heat  resistance  of  the  bulkhead,  and 
lightening  the  bulkhead.  The  bulkhead  structure  and  function,  which  are  the 
most  important  technologies,  were  studied  in  this  research  work.  A  porous 
system,  a  90°  opening  and  rotating  door  system,  a  slide  opening  system,  a  heat 
insulating  material  peeling  system,  etc.,  were  studied  in  a  pre-experiment. 
We  have  come  to  the  conclusion  that  current  design  methods  and  materials  are 
insufficient,  except  for  the  porous  system,  in  terms  of  reliability  of 
operation,  simplicity  and  weight.  As  a  result,  we  have  determined  to  adopt 
the  porous  system.  We  have  made  a  bulkhead  testpiece  suitable  in  size  for  an 
experimental  scale,  have  incorporated  this  bulkhead  testpiece  in  a  wall 
thickness-type  small-sized  rocket  motor  with  a  diameter  of  80  mm,  and  have 
conducted  an  experiment  on  the  combustion  of  this  rocket  motor.  Figure  7 
shows  a  structural  drawing  of  the  rocket  motor  used  in  this  experiment. 


1  Pickup  for  measurement  of  second  stage  pressure 

2  Second  stage  igniter 

3  Second  stage  propellant 

4  Bulkhead 

5  First  stage  propellant 

6  First  stage  igniter 

7  pressure  leading  hole  for  measurement  of  the 
first  stage  pressure 

8  Nozzle 

Figure  7.  Experimental  Device 


Figure  8.  Structural  Drawing  of  Experimental  Bulkhead 
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Figure  8  shows  a  structural  drawing  of  the  bulkhead  used  in  the  experiment. 
This  bulkhead  is  a  stainless  disk  with  a  diameter  of  80  mm.  Twenty  holes  with 
a  diameter  of  4.5  mm  are  made  in  the  bulkhead;  the  total  area  of  the  holes  is 
2.8  times  that  of  the  nozzle  straw  area.  Therefore  the  speed  when  these  holes 
in  the  bulkhead  are  open  and  the  combustion  gas  passes  through  them  will  be 
low-subsonic,  about  Mach  0.2,  and  will  not  result  in  any  pressure  loss  that 
would  lower  the  performance  of  the  rocket  motor.  Ceramic  balls  with  a  diameter 
of  6  mm  were  fastened  to  the  first  stage  combustion  chamber  side  of  these 
holes  in  order  to  prevent  the  heat  and  combustion  gas  generated  during  first 
stage  combustion  from  penetrating  the  second  chamber.  This  diameter  is  about 
half  that  of  the  nozzle.  These  ceramic  balls  are  designed  so  that  when  the 
second  stage  propellant  is  ignited,  the  gas  pressure  will  separate  them  from 
the  bulkhead  and  they  will  be  discharged  to  the  outside  of  the  nozzle  through 
the  first  stage  combustion  chamber,  which  is  already  empty. 

Figure  9  [not  reproduced]  is  a  photograph  of  the  bulkhead  used  in  this 
experiment.  It  was  taken  from  the  side  of  the  first  stage  combustion  chamber. 
The  white  sections  are  the  ceramic  balls  affixed  to  the  first  stage  combustion 
chamber  side.  A  mixture  of  epoxy  resin  and  titanium  oxide  was  applied  as  a 
heat  resistant  material  around  the  ceramic  balls  in  order  to  raise  the  heat 
resistance  of  the  bulkhead. 

Table  1  and  Figure  10  show  the  composition  and  shape  of  the  propellant  used  in 
the  experiment,  respectively.  The  internal  both-end-face  combustion  system  was 
adopted  as  the  propellant  combustion  system  to  keep  combustion  chamber  pres¬ 
sure  as  constant  as  possible.  The  theoretical  flame  temperature  was  2,900  K. 

Table  1.  Composition  of  Sample  Propellant 

Kind  of  propellant  R-1 


Component  Particle  diameter  (/im)  (Percent) 


Figure  10.  Shape  of  Sample  Propellant 


Carbon  graphite  was  used  as  the  nozzle  material.  The  nozzle  throat  diameter 
(Dt)  was  13.58  mm  and  the  opening  ratio  (e)  was  2.0. 


An  igniter  installed  in  a  hole  in  the  propellant  was  used  to  activate  the 
rocket  motor.  An  ignition  code  for  the  first  stage  combustion  chamber  and 
another  for  the  second  stage  combustion  chamber  were  passed  from  the  nozzle 
and  the  forward  cap,  respectively.  The  HTPB/AP  composite  propellant  was  cut 
into  5  mm  cubes  and  these  small  cubes  were  mixed  with  6  g  of  black  powder. 
This  mixture  was  used  as  a  powder  for  the  igniter. 

The  measurement  variable  was  determined  to  be  the  pressure  in  each  combustion 
chamber.  The  first  stage  and  the  second  stage  combustion  chamber  pressures 
were  measured  with  a  strain  gauge  pressure  pickup  via  a  pressure  leading  hole 
in  the  nozzle  cap  and  a  forward  cap,  respectively. 

4.  Results  and  Discussion 

Figure  11  shows  the  relationship  between  combustion  time  and  combustion 
chamber  pressure  obtained  from  the  experiment. 


First  stage  combustion  time:  5  seconds 
Combustion  interruption  time: 13  seconds 
Second  stage  combustion  time:  3  seconds 


I - _J - 1 _ I _ I _ I _ I 

0  2  4  6  1  8  2  0  2  2 


Time  (S6C) 


Figure  11.  Relationship  Between  Combustion  Pressure  and  Time 


First,  the  first  stage  igniter  was  activated  to  burn  the  first  stage 
propellant.  This  propellant  was  burned  for  5  seconds  at  20  kgf/cm^  without  any 
ignition  delay. 


It  can  be  determined  that  the  ceramic  balls  adhered  to  the  bulkhead 
effectively,  i.e.,  they  completely  shut  off  the  first  stage  combustion  gas  and 
heat.  This  can  be  verified  by  the  fact  that  the  pressure  in  the  second  stage 
combustion  chamber  did  not  increase  during  the  5-second  burn. 

After  the  second  stage  combustion  was  completed,  and  after  13  seconds  passed, 
the  second  stage  igniter  was  activated.  As  soon  as  the  second  stage  propellant 
was  ignited,  the  ceramic  balls  were  separated  from  the  bulkhead  by  the 
combustion  gas  pressure  and  were  discharged  from  the  nozzle  via  the  flow  of 
combustion  gas.  There  was  a  possibility  of  combustion  chamber  pressure  being 
generated  from  the  discharge  of  the  ceramic  balls ,  but  no  rise  or  peak  of 
harmful  pressure  like  the  combustion  chamber  pressure  was  detected. 

The  second  stage  propellant  was  burned  for  3  seconds  at  40  kgf/cm^  without  any 
ignition  delay.  There  was  no  difference  between  the  pressure  of  the  first 
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stage  combustion  chamber  and  that  of  the  second.  Therefore,  it  can  be 
concluded  that  the  combustion  gas  flowed  from  the  second  stage  combustion 
chamber  to  the  first  without  any  pressure  loss  caused  by  the  bulkhead.  This 
experiment  indicates  that  there  was  no  interference  (influence  by  heat 
transfer,  etc.)  of  the  first  stage  combustion  chamber  pressure  to  the  second 
during  the  interruption  of  combustion  for  at  least  13  seconds,  and  that  the 
combustion  interruption  time  can  be  selected  arbitrarily  for  an  interval  of  13 
seconds  or  less. 

After  the  operation  of  the  rocket  motor  was  completed,  it  was  disassembled  for 
operation.  Figure  12  [not  reproduced]  shows  the  condition  of  the  bulkhead. 
All  the  ceramic  balls  were  separated  from  the  bulkhead,  the  degree  of  burn 
damage  to  the  bulkhead  was  very  low,  and  the  holes  on  the  second  stage 
combustion  chamber  side  were  burned  and  slightly  damaged.  There  were  no 
ceramic  balls  in  the  combustion  chamber;  all  the  ceramic  balls  were  discharged 
through  the  nozzle  to  the  outside.  Also,  neither  breakage  nor  damage  was  seen 
on  the  nozzle.  Therefore,  it  can  be  concluded  that  all  the  ceramic  balls 
passed  smoothly  through  the  nozzle. 

Based  on  the  above,  it  can  be  said  that  the  rocket  motor  used  in  this 
experiment  achieved  thrust  control  without  lowering  combustion  performance. 

5.  Conclusion 

(1)  The  bulkhead  breaking  system  used  in  this  research  study  had  no  influence 
on  the  combustion  pressure  and  temperature  in  the  second  stage  combustion 
chamber  during  first  stage  combustion.  There  was  no  damage  to  the  nozzle, 
etc.,  due  to  fragments  from  the  bulkhead,  nor  was  there  any  decrease  in  the 
performance  of  the  rocket  motor. 

(2)  No  pressure  loss  was  seen  at  the  bulkhead  through  which  the  combustion 
gas  of  the  second  stage  propellant  passed. 

(3)  Stable  combustion  was  achieved  by  a  thrust-control  rocket  motor  using  a 
bulkhead  based  on  this  bulkhead  breaking  system.  Also,  we  were  able  to 
control  the  operation  of  this  rocket  motor  for  a  specified  time  interval 
(about  10  seconds,  virtually  the  assumed  time). 

6.  Future  Research  Subjects 

In  the  future,  it  will  be  necessary  to  carry  out  research  in  the  following 
areas  by  using  an  actual  missile:  1)  lightening  the  bulkhead;  2)  testing  and 
evaluating  heat  resistance;  and  3)  the  ignition  system  for  each  stage. 
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